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ABSTRACT 


The  phenomenon  of  disturbances  caused  by  electric 
arc  discharge  in  a  supersonic  flow  was  studied  both  analyti¬ 
cally  and  experimentally.  The  purpose  was  to  study  the 
feasibility  of  using  this  concept  for  attitude  control  of 
space  vehicles.  A  possible  application  for  rocket  propelled 
vehicles  was  selected  for  the  test  condition  of  an  experi¬ 
mental  program.  Fair  agreement  was  obtained  between  the 
test  results  and  the  analytical  solutions.  It  was  found  that: 

1.  A  local  static  pressure  rise  can  be  created  by 
the  electric  discharge  in  a  supersonic  stream. 

2.  The  electric  field  to  initiate  discharge  was 
about  2  x  103  volts/in  for  nitrogen  stream  at 
Mach  number  2.5  with  static  pressure  of  2  psia 
and  stagnation  temperature  of  4000°R. 

3.  The  arc  discharge,  without  a  magnetic  field,  in 
a  supersonic  flow  of  Mach  number  2.5  can  be 
maintained  and  is  stable. 

4.  The  affected  area  is  small  which  resulted  in 
small  side  force. 

5.  The  generated  side  force  per  unit  electric  power 
input  is  small  in  view  of  the  available  electric 
power  supply. 

The  use  of  an  electric  discharge  for  practical  thrust 
vector  control  of  rocket  motors  is  only  feasible  if  the  power 
supply  can  be  improved.  This  concept,  however,  proved  to  be 
physically  sound  and  may  be  very  important  for  applications 
in  nuclear  rockets  where  electrical  energy  is  already  avail¬ 
able  on  board,  or  to  the  space  mission  re-entry  vehicles 
v/hich  the  needed  power  has  been  generated  and  stored  before 
re-entry . 
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Attitude  control  is  needed  for  both  rocket  propelled  and 
re-entry  vehicles.  The  first  device  to  achieve  this  purpose 
was  to  control  the  jet  var.es  or  to  gimble  the  rocket  exhaust 
nozzles  by  mechanical  means.  The  scheme  ci  secondary  fluid 
injection  was  introduced  in  1952  (Hef .  1)  by  using  mass  addition 
into  the  rocket  exhaust  nozzle  to  create  the  desired  side  force 
for  thrust  vector  control.  Analytical  studies  made  in  early 
1950 's  (Reis.  2,  3,  and  4)  provided  several  different  approaches 
for  predicting  the  generated  side  force,  and  also  revealed  the 
possible  applications  for  re-entry  vehicles,  Kaufman  (Hef.  5) 
in  1957  shoved  his  experimental  results  in  hypersonic  regime. 

The  equivalence  principle  between  heat  and  mass 
sources  in  high  speed  fleas  was  noticed  by  Sicks  (Ref.  6)  in 
1950  and  by  Chu  (Ref.  7)  in  1955,  and  further  developed  by 
Willmarth  (Ref.  S)  in  1957.  This  concept  leads  to  the  present 
investigation  of  using  electric  arc  discharge  to  produce  the 
heat  source  for  vehicle  attitude  control. 

Compared  with  the  existing  techniques,  the  electric 
discharge  concept  premises  the  following  inherent  advantages: 

i)  It  has  no  moving  parts, 

ii)  Plumbing  system  is  not  required, 

iii)  Instantaneous  response  rate  appears  possible, 

iv)  Since  the  system  is  electrical,  it  may  be  more 

compatible  with  the  guidance  and  control  systems. 

The  possible  disadvantages  may  be  as  f  ol  lews: 

i)  It  may  result  in  a  heavier  system  compared  to  that 
of  the  secondary  injection  scheme, 
ii)  There  may  be  material  problems  associated  with  high 
arc  temperatures. 

This  report  is  a  study  of  the  feasibility  of  using  the 
new  method  to  steer  a  vehicle  during  its  flight.  The  proposed 
technique  is  an  attempt  to  obtain  a  positive  and  instantaneous 
control  of  the  system  by  means  of  discharging  electric  energy 
locally  into  the  supersonic  stream  as  shown  in  Fig.  1. 

For  rocket  applications,  the  electric  discharge  is 
accomplished  by  an  arc  located  at  one  side  of  the  nozzle  wall, 
which  causes  a  static  pressure  rise  in  the  vicinity  of  the 


discharge.  This  asymmetric  pressure  rise  inside  the  nozzle 
exit  cone  results  in  a  side  force  exerted  against  the  nozzle 
wall.  A  similar  phenomenon  has  been  observed  in  the  secondary 
fluid  injection  scheme  for  thrust  vector  control  of  rockets. 

The  same  concept  may  also  be  applicable  to  the  re-entry  vehicles. 

In  order  to  gain  some  insight  o':  the  phenomenological  aspects 
of  this  concept,  a  feasibility  study  has  been  conducted  both 
theoretically  and  experimentally.  The  magnitude  of  the  side 
force  was  predicted  by  three  different  analytical  models  as 
discussed  in  Section  II.  The  electric  discharge  and  its  cor¬ 
responding  pressure  rises  were  investigated  by  a  limited 
experimental  program  as  discussed  in  Section  III. 
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SECTION  II 


ANALYTICAL  INVESTIGATION 


2.1  THEORETICAL  MODELS 

The  phenomena  associated  with  electric  discharge  in  super¬ 
sonic  flow  (Refs.  9,  10,  and  11)  are  very  complex.  The  flow 
field  is  affected  by  the  electric  discharge,  which  in  turn  will 
be  influenced  by  the  flow.  The  side  force,  produced  by  the 
interaction  of  the  flow  field  and  the  electric  discharge,  can 
only  be  analyzed  based  on  simplified  models  at  the  present 
tine.  The  validity  of  these  simplifications,  however,  will  be 
verified  by  experimental  investigation. 

From  experience  in  secondary  fluid  injection  into  a  super¬ 
sonic  flow,  it  was  found  that  both  the  ,fBlast  Wave  Analogy," 
(Hef.  4),  and  the  "Linearized  Snail  Disturbance  Theory"  (Ref.  2) 
gave  a  fairly  good  gross  agreement  with  the  measured  side  force, 
la  addition  to  the  above  two  models,  a  "Plane  Flame  Analogy" 
was  also  employed  to  predict  the  magnitude  of  the  generated 
side  force.  All  three  models  gave  the  same  order  of  magnitude 
results  and  are  discussed  as  follows. 

2.1.1  Blast  gave  Analogy 

The  flow  field  induced  by  a  blast  wave  in  a  quiescent 
atmosphere  was  studied  extensively  by  various  authors  (Beife.  12, 
13,  14,  and  15).  Lees  (Ref.  16)  observed  that  the  bow  shock 
wave  formation  due  to  blunt  nosed  bodies  in  a  hypersonic  flow 
is  equivalent  to  a  constant  energy  supply  in  the  form  of  a 
body  drag.  He  successfully  related  the  unsteady  phenomenon  of 
a  sudden  energy  release  into  a  quiescent  flow  to  the  steady 
phenomenon  of  an  aerodynamic  body  in  hypersonic  flow.  The 
electric  arc  discharge  in  a  supersonic  flow  can  be  considered 
as  a  high  thermal  energy  release  into  a  high  velocity  gas 
stream. 

The  combined  effects  of  both  the  thermal  energy  and  the 
free  stream  velocity  makes  it  possible  to  apply  the  concept 
of  the  blast  wave  analogy  for  analyzing  electric  discharge 
and  flow  field  interactions.  It  is  obvious  that  the  accuracy 
obtained  from  this  analogy  can  be  improved  as  the  thermal  energy 
or  the  free  stream  velocity  increases.  Roman  (Ref.  17)  also 
observed  that  an  electric  arc  confined  by  a  magnetic  field 
behaves  like  a  heated  solid  body  in  a  flow  field.  This  observa¬ 
tion  lends  a  strong  support  of  the  validity  of  the  blast  wave 


analogy  applying  to  electric  arc  discharge  into  a  supersonic 
flow  field. 

In  order  to  see  the  nature  of  the  problem  and  to  determine 
the  equivalent  energy  source  strength,  it  is  assumed  that  a 
large  amount  of  heat  is  released  on  a  flat  plate  which  has  a 
zero  angle  of  attack  in  a  uniform  supersonic  stream.  The  energy 
dissipation  in  the  viscous  boundary  layer  is  small  compared  to 
the  energy  release  due  to  the  blast  and  was  ignored.  The  net 
result  of  interaction  between  the  flow  field  and  the  released 
energy  can  then  be  adequately  treated  by  the  blast  wave  theory. 

If  the  energy  is  released  at  a  point  at  a  constant  rate, 
the  pressure  wave  will  be  spherically  symmetric.  This  corresponds 
to  a  point-blast  case.  A  uniform  flow  of  velocity  u  super¬ 
imposed  on  the  point  energy  source  will  result  in  an  axially 
symmetrical  steady  flow.  In  other  words,  the  phenomenon  is 
equivalent  to  that  of  the  point  energy  source  moving  with  a 
speed  of  u  in  the  upstream  direction  with  respect  to  a  quies¬ 
cent  flow.00  This  constitutes  the  strength  of  the  line  energy 
source  per  unit  length,  E,  and  can  be  expressed  as 
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CO 


(1) 


where  P.  is  the  total  energy  input  per  unit  time,  or  the  pov/er 
input.  1The  above  relation  may  be  understood  in  a  different 
manner;  the  total  power  input  to  the  flow  has  to  be  equal  to 
the  strength  of  a  line  source  carried  with  the  characteristic 
velocity  in  question,  namely,  the  free  stream  velocity  u^. 

The  above  equation  is  derived  without  considering  the  wall 
effect.  If  the  source  is  on  the  wall,  and  without  considering 
the  boundary  layer,  this  means  that  only  the  half-space  is 
affected  by  the  above  energy  release;  that  is,  half  of  the 
strength  will  produce  the  same  net  result  in  the  half-space. 
Hence,  in  the  flat  plate  case,  the  energy  per  unit  length  is 
given  by 
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The  first  order  similarity  solution,  from  blast  wave  theory, 
for  a  strong  shock  wave  is  given  by  Sakurai  (Ref.  7.3)  and 
Sedov  (Ref.  14)  for  the  effected  area  as 


(—  )2(— )2  =  i  (7)  (3) 
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and  for  the  magnitude  of  pressure  as, 
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where  R  is  defined  as, 
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and  a  is  the  free  stream  speed  of  sound,  R  is  the  shock  radius 
u  =  dS/dt  is  the  shock  speed,  jQ  is  a  constant  depending  on  the 
ratio  of  specific  heats  7,  p  is°the  pressure,  t]  =  r/R,  with  r 
the  distance  from  the  center  of  explosion,  and  g(T})  is  a  non- 
dimensional  pressure  given  as: 
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If  the  energy  is  released  adjacent  to  a  wall,  the  force 
exerted  on  the  wall  can  be  obtained  by  integrating  the  pressure 
rises  over  the  effective  area  as  (Fig.  2), 


LR 

F  =  2  //  (p  -  p)drdx  (7) 
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where  L  is  an  arbitrary  length  toward  downstream,  and  R  is  to 
the  shock  position.  Broadwell  (Ref.  4)  successfully  showed 
that  Equation  (7)  can  be  integrated  for  the  maximum  side  force 
by  employing  the  relations  of  Equations  (3),  (4),  (5),  and  (6) , 
resulted  in: 


number  and  o(y)  is  a  constant 


where  is  free  stream  Mach 
depending  on  the  value  of  7, 


o(y)  =0.10 
o(y)  =  0.14 
o(y)  =0.17 
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for 

7  =  1.4 

The  side  force  F  exerted  on  the  plate  due  to  a  certain 
power  input  P^n,  in  tie  present  problem,  then  becomes 
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or,  the  side  force  per  unit  power  input  is 
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In  terns  of  the  line  energy  strength, 
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Equation  (10c)  indicates  that  the  side  force  generated  is 
linearly  proportional  to  the  free  stream  Mach  number  for  the 
same  line  energy  strength.  Hence,  it  will  be  more  beneficial 
to  have  the  energy  released  to  a  higher  Mach  number  stream 
than  the  lower  Mach  number  supersonic  stream.  A  plot  of  the 
generated  side  force  vs.  free  stream  Mach  number  is  shown  in 
Fig.  3. 


2.1.2  Supersonic  Small  Disturbance  Theory 


The  present  problem  of  heat  release  to  a  supersonic  flow, 
if  viewed  from  the  other  angle,  namely,  the  entire  flow  field, 
may  be  said  to  only  be  a  local  phenomenon.  The  energy  addition 
only  constitutes  a  perturbation  to  the  main  flow,  hence 
linearized  theory  can  be  used. 

Linearized  supei’sonic  small  disturbance  theory  has  been 
employed  in  Wu,  Chapkis  and  Hager's  (Ref.  2)  analysis  and  showed 
a  good  agreement  in  predicting  the  order  of  magnitude  of  side 
force  generation  by  fluid  injection.  Broadwell  (Ref.  4)  also 
pointed  out  that  the  linearized  theory  should  give  the  right 
order  of  magnitude  of  the  side  force.  However,  Broadwell  added 
one  additional  term  to  the  side  force  equation  proposed  by 
Wu,  et  al .  He  argued  that  the  supersonic  small  disturbance 
theory  implicitly  included  the  mass  source  entered  to  the  free 
stream  with  an  x-component  velocity  of  the  same  magnitude  of 
the  free  stream.  The  addition  of  that  term  made  the  result  of 
the  equation  by  an  order  of  magnitude  greater  than  the  original 
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one-  The  disturbance  created  by  the  plane  source  flow  in  a 
supersonic  stream  is  caused  only  by  the  "interaction"  of  the 
main  and  the  source  streams.  There  is  no  need  to  add  another 
term  to  take  care  of  the  artificial  x-component  of  ihe  velocity 
since  the  disturbance  created  by  the  plane  source  flow  is  a 
net  interaction  phenomenon.  A  detailed  discussion  of  Sroadwell’s 
additional  term  is  included  in  Appendix  I . 

The  analysis  of  a  plane  mass  injection  into  a  supersonic 
stream  on  a  flat  plate  was  performed  by  Wu,  et  al.  The  pressure 
coefficient  C  due  to  the  mass  injection  disturbance  is 
analogous  to  the  disturbance  created  by  the  half -wedge  angle 
6  (Fig.  4) 
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where  vu  is  the  vertical  disturbance  velocity.  If  the  secondary 
mass  flow  rate  m^,  over  the  area  A,  is  small  as  compared  to 
the  mainstream  mass  flow  rate,  then  the  vertical  disturbance 
velocity  can  be  related  to  the  injected  mass  flow  rate-  as 
follows: 
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Hence,  the  pressure  disturbance  by  the  mass  injection  can  be 
written  as 
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The  side  force  exerted  on  tbe  flat  plate  is  therefore 
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Assuming  that  the  disturbances  were  created  by  a  mass  source, 
a  heat  source  ard  a  body  force,  Willmarth  (Ref.  8)  derived  the 
governing  flow  equation  of  the  x-comoonent  as  follows  (ApDendix 
ID, 
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where  n  is  the  mass  source  per  unit  volume,  q  is  the  heat  source 
per  unit  volume,  and  f  .  is  the  x-component  of  the  body  force  per 
unit  volume.  It  can  be  seen  that  the  mass  source  and  heat 
source  can  be  related  as 


m 


(16) 


where  m  and  Q  are  the  total  quantities  of  the  rate  of  mass  and 
heat  additions  respectively.  This  analogy  was  also  observed  by 
Hicks  (Ref.  6)  and  Chu  (Ref.  7). 

Employing  Equation  (16) ,  the  generated  side  force  on  the 
wall  as  expressed  in  Equation  (14)  can  be  written  as 
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the  side  force  per  unit  power  input  can  be 


Since  Q  =  P^, 
expressed  as 
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It  is  surprising  to  find  out  that  the  asymptotic  fora  of 
the  above  equation  for  large  M_  is  similar  to  that  of  the  blast 
wave  theory,  although  linearized  assumption  should  not  apply 
to  high  Mach  number  flows.  Nevertheless,  both  theories  give 
the  relation  that 
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The  plot  of  generated  side  force  vs.  free  stream  Mach  number 
is  shown  in  Fig.  5. 

2.1.3  Plane  Flame  Analogy 


There  is  another  possible  way  to  calculate  the  magnitude 
of  the  side  force  by  a  simple  plane  flame  analogy.  It  is  assumed 
that  the  heat  is  added  into  an  inviscid  supersonic  stream  through 
a  plane  thin  flame  as  shown  in  Fig.  6.  This  approach  is  similar 
to  Willmarth's  (Ref.  8)  flame  method  except  that  the  solutions 
of  the  flow  parameters  can  be  uniquely  determined. 

The  analysis  of  the  flow  field  resulting  from  the  heat 
addition  through  a  plane  flame  can  be  greatly  simplified  by 
making  the  following  assumptions: 

i)  The  gas  is  thermally  and  calorically  perfect, 

ii)  Gas  properties  do  not  change  across  the  flame, 

iii)  Plane  flame  has  a  negligible  thickness,  standing 
at  an  angle  a  to  the  plate  (Fig.  6). 

iv)  Heat  transfer  occurs  only  when  the  flow  passes 
this  flame. 


Across  a  flame  sheet,  the  normal  ccmpememt  ci  the  flew 
velocity  ma. y  either  increase  cr  decrease  depending  ca  whether 
tfce  flare  is  undergoing  a  deflagration  cr  detonation  process. 

In  either  case ,  the  velocity  ccnrpcneni  parallel  to  the  flans 
is  unaffected.  Therefore,  a  deflagration  flame  will  torn  the 
flew  toward  the  plate,  while  a  detonation  flame  will  turn  the 
flew  away  from  the  plate.  However,  the  boundary  condition 
requires  that  the  flew  must  be  parallel  to  the  plate,  therefore 
the  presence  of  a  second  wave  is  necessary.  There  are  two 
possible  waves  that  nay  be  present,  sanely ,  a  shock  wave  which 
always  decelerates  the  coma!  (to  shock  wave)  component  of 
the  velocity,  hen 'e  turns  the  flew  away  from  the  plane,  cr 
expansion  waves  which  always  turns  the  flew  coward  the  plate, 
in  the  configuration  as  shewn  in  Fig.  S.  Therefore,  we  may 
have  eight  possible  flew  patterns,  as  listed  in  Table  I. 


There  are  four  cases  where  the  boundary  conditions  are 
satisfied.  Case  (1)  (a  shock  wave  followed  by  a  deflagration 
flame)  does  not  violate  any  basic  physical  or  mathematical 
laws,  therefore  its  existance  is  possible.  Case  (4)  (an 
expansion  wave  followed  by  a  detonation  flame)  is  ruled  cut 
because  of  the  following  consideration.  In  supersonic  flow 
(here  only  the  velocity  component  normal  to  these  waves  are 
considered,  the  parallel  component  may  be  eliminated  if  the 
observer  moves  with  that  speed),  all  disturbances  propagate 
downstream  only.  Tbe  existence  of  an  expansion  wave  upstream 
Of  the  detonation  fl2me  requires  the  signals  propagate  ap>- 
stream.  Therefore,  Case  (4)  can  not  occur.  Case  (5)  Is  a 
deflagration  flame  followed  by  a  shock.  However,  there  is  no 
deflagration  flame  in  supersonic  flow,  this  rules  out  the 
possibility  of  Case  (5).  Finally,  we  shall  consider  Case  (8) 
(detonation  followed  by  an  expansion) .  Expansion  waves  can 
only  occur  in  supersonic  flow.  In  the  present  case  it  re¬ 
quires  the  flow  be  supersonic  (the  normal  component)  after 
the  detonation  flame.  This  type  of  flame  Is  called  weak 
detonation  which  h2s  not  been  observed  In  reality.  Therefore, 
Case  (8)  does  not  exist.  The  only  condition,  which  both 
satisfies  the  boundary  condition  and  is  physically  possible, 
is  a  shock  wave  followed  by  a  deflagration  flame  (Case  (1)). 
Therefore,  the  solution  is  uniquely  determined. 
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The  continuity,  mcnentun  and  energy  equations,  together 
with  the  eqnatioa  of  state,  for  the  flow  normal  to  the  shcck 
wave  are: 


and 


c  u  sin  p  =  £  u  sin(p  -  6) 
12.  2  2 
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Those  equations  for  the  flow  normal  to  the  flame  are: 


p  u  sin(a  -  6)  =  p  u  sin  a  (23) 

2  2  3  3 


p  +  p  u2  sin2(a  -  5)  =  p  +  p  u2sin2a  (24) 

2  2  2  3  3  3 


C  T  +  -  u2  +  AQ  =  C  T  +  --  u2  (25) 
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(26) 


p  /p  T  =  p  /p  T 

2  2  2  3  3  S 


Trie  Chapnan-Jcruguet  deflagration  process  across  the  flare  is: 


u  sin  a  =  ("vST  )2 
3  3 


(27) 


In  addition,  two  equations  expressing  the  fact  that  the  tan¬ 
gential  velocity  does  not  change  across  the  shock  and  the  flare 
nay  be  written  as: 


u  cos  0  =  u  cos  O  -  5)  (28) 

1  2 


u  cos  (a  -  6)  =  u  sin  a  (29) 

2  3 


where  p  is  the  density,  u  the  velocity,  T  the  temperature,  p 
the  pressure,  AQ  the  heat  addition  per  unit  mass  of  gas,  {3 
the  oblique  shock  wave  angle,  a  the  flame  angle,  6  the 
deflection  angle  across  the  shock  wave  and  the  subscripts, 
x  2  3  are  indicating  the  regions  respectively,  as  shown  in 
Fig-  6.  Now  we  have  eleven  equations  with  eleven  unknowns. 

This  problem  is  solveable.  However,  due  to  the  complexity  of 
algebra  involved  in  simplifying  the  above  equations ,  it  is  felt 
that  the  best  method  is  to  solve  numerically  and  iterate  the 
results  with  the  help  of  tabulated  numerical  tables  for  the 
static  temperatures  and  pressures  cf  one-dimensional  heat 
addition. 

The  generated  side  force  per  unit  power  input  can  be  cal¬ 
culated  by  considering  the  geometrical  relationships  as  shown 
in  Fig.  6.  The  total  power  input  through  the  flame  is  p  u  AQ, 
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hence  the  side  force  per  unit  pover  input  is 
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The  side  force  per  unit  pover  input  calculated  according 
to  Equation  (30)  is  in  agreement  vith  the  order  of  magnitude  to 
that  of  the  blast  vave  theory.  Side  forces  evaluated  by  both 
the  blast  wave  theory  and  the  flare  theory  appeared  to  be  less 
than  that  of  the  linearized  theory.  This  is  to  be  expected 
because  the  linearized  theory  implicitly  assumed  reversible 
heat  addition  and  should  give  the  minimr-)  loss.  Hovever ,  in 
reality,  the  blast  vave  theory  and  the  flame  theory  should 
give  better  results  due  to  the  noa— linear  effect  caused  by 
the  large  amount  of  heat  addition. 

2.1.4  Summary  of  Theoretical  Studies 

Frcn  the  above-mentioned  three  approaches ,  one  can  Dredict 
the  magnitude  of  the  side  force  generated  by  heat  addition  from 
a  side  vail  to  the  main  stream.  For  a  special  case  vhich  is 
investigated  later  experimentally,  the  folloving  conditions  are 
used: 
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the  specific  side  force  per  unit  power  input  F  /P.  for  the 
three  approaches  are  as  follows:  s  xa 


. . 

Blast  Wave 

Linearized 

Plane  Flame 
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<\.  .035 

where  F  /P.  has  the  unit  of  lbs/kw,  or  approximately,  Ibs/BTU/sec. 
All  of  fhe1three  different  approaches  predict  the  side  forces  of 
the  same  order  of  magnitude.  Blast  wave  and  linearized  theories 
give  the  solution  explicitly,  while  the  plane  flame  theory  needs 
iteration.  As  expected,  the  blast  wave  theory  gives  the  better 
agreement  than  the  linearized  theory  as  compared  to  our  experi¬ 
mental  results  which  will  be  presented  in  Section  III. 

The  linearized  theory  predicts  the  highest  specific  side 
force  per  unit  power  input.  This  is  so  because  the  small  dis¬ 
turbance  theory  represents  the  condition  of  minimum  loss.  There¬ 
fore,  the  value  may  be  considered  as  the  upper  limit  of  side 
force  obtainable  by  the  arc  discharge  in  a  low  supersonic  flow. 

For  a  high  supersonic  flow  or  large  energy  addition,  the  linearized 
theory  breaks  down  due  to  its  inherent  assumption  which  we  made 
in  the  derivation  of  the  governing  equations. 

The  blast  wave  theory,  on  the  other  hand,  can  not  extrap¬ 
olate  to  the  low  Mach  number  flow  or  low  energy  addition  case. 

This  is  due  to  the  assumption  of  a.  strong  shock  wave  boundary 
condition  employed  in  its  original  derivation. 


2.1.5  Comparison  of  Heat  Addition  and  Mass  Injection  in 
Generating  Side  Fcrce 


From  the  above  discussions,  it  is  clear  that  both  mass 
injection  and  electrical  energy  addition  are  capable  to  generate 
side  force  locally  in  a  supersonic  stream.  Purely  from  a  the¬ 
oretical  point  of  view,  it  is  difficult  to  determine  which  method 
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is  superior.  Depending1,  on  the  quantities  of  injection,  both 
schemes  may  yield  comparable  magnitudes  of  side  force.  There¬ 
fore,  the  selection  of  either  method  for  a  specific  mission, 
depends  very  much  on  the  conditions  such  as  the  state  of  the 
art  on  power  source,  the  material  capability,  the  hardware 
designs,  etc.  In  some  cases,  there  is  electrical  power  avail¬ 
able  on  board  such  as  nuclear  propelled  vehicles,  the  electrical 
energy  addition  through  the  arc  discharge  technique  appears  to 
be  superior,  because  there  is  no  need  to  carry  the  redundant 
auxiliary  tank  for  the  mass  injection.  On  the  other  hand,  the 
arc  discharge  method  depends  on  the  specific  energy  of  the 
power  source  if  tb/e  power  source  needs  to  be  carried  on  board. 

The  specific  energy  of  the  power  source  is  defined  as  the  ratio 
of  the  energy  of  the  source  to  its  weight . 

From  the  airborne  weight  penalty  point  of  view,  the 
performance  of  arc  discharge  may  be  best  analyzed  by  using  the 
parameter  of  side  force  specific  impulses,  Isps*  This  parameter 
is  defined  as  the  side  force  generated  per  unit  weight  consumption. 

For  the  mass  injection  case, 
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sps^m 


(31) 


and  for  the  electric  discharge  case, 
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where  w  =  flow  rate  of  secondary  injectant,  W  =  battery,  or 
electric  power  sour: e ,  consumption  rate  definld  as 
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and  E  =  battery  specific  energy,  usually  is  given  by 
watt-lP/lb.  Hence, 


(Isps>e 


F  E 
s  sp 


£"in 


(34) 


The  side  force  specific  impulse  can  be  evaluated  by  the  blast 
wave  theory  and  the  linearized  theory,  respectively  * 

2 .1.5.1  Blast  Wave  Theory 

The  side  force  generated  by  the  thermal  energy  release  has 
been  derived  as 


Fseb  * 


Pin 


u 


(35) 


Hence,  the  specific  impulse  may  be  written  as 


(I 


sps^eb 


=  d(7)Moo  -2£  (36a) 

00  u 

a? 


by  employing  the  conventional  unit,  i.e.,  E  in  watt-hr/lb,  u 

in  ft/sec,  I  in  sec,  the  above  equation  Becomes 
sps 


(Isps>eb  -  2655 


E 

SP 

u 


(36b) 
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2.1 >5. 2  Linearized  Snail  Disturbance  Theory 

The  same  argument  may  be  also  applied  to  the  linearized 
theory.  The  expression  for  specific  impulse  of  electric  dis¬ 
charge  based  on  linearized  theory  is  as  follows: 


^sps^et 


2655  ~ 


M2 

03 


■V»£  -  1 


sp 


u 


(37) 


The  theoretical  performance  of  employing  electric  discharge 
in  generating  side  force  by  the  above  two  theories  are  shown  in 
Fig.  7.  The  side  force  specific  impulse  is  linearly  related  to 
the  specific  energy  of  the  power  source.  At  the  present  state, 
the  specific  impulse  of  using  electric  discharge,  by  carrying 
chemical  batteries  on  board,  is  significantly  lower  than  that 
of  the  secondary  mass  injection  because  of  the  low  specific 
energy  values.  However,  the  electric  discharge  technique  will 
appear  to  be  superior  than  the  secondary  fluid  injection  scheme, 
if 

i)  the  electrical  source  is  available  on  board,  such  as 
nuclear  rockets, 

ii)  storage  of  solar  energy  is  possible  during  the  long 
time  space  mission  vehicles  in  its  re-entry  phase,2 

iii)  further  improvement  and  development  of  chemical 
batteries  are  possible,3  and 

iv)  development  of  new  type  power  source  like  MHD  generators4 
becomes  available. 


2At  the  present  capability,  a  round  trip  to  outer  planets 
takes  several  years  in  space.  In  addition,  according  to  Cheng 
(Ref.  21)  the  arc  discharge  applied  to  a  high  Mach  number  low 
density  flow  may  result  in  some  other  effects  which  will  help 
in  generating  side  force. 

3It  is  felt,  with  the  present  technology  on  chemical  batteries, 
that  the  limit  of  specific  energy  required  from  this  system  is 
below  that  required  to  make  thrust  vector  control  by  electric 
arc  discharge  competitive. 


4Airborne  type  MHD  generators  are  in  the  development  stage  at  the 
J.  B.  Dicks  and  Associates,  Tullahoma,  Tennessee  (see  detail  in 
Appendix  III). 
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2.2  ELECTRIC  ARC  CHARACTERISTICS 


An  analytical  method  is  not  available  for  predicting  either 
the  breakdown  voltage  or  the  operating  characteristic  of  an  elec¬ 
tric  arc  with  sufficient  confidence.  For  the  configurations 
that  the  electrodes  are  mounted  flush  on  the  same  surface  with- 
oxrt  facing  at  each  other,  there  is  no  information  available 
even  to  serve  as  a  guide  for  an  intelligent  guess.  It  is  felt 
that  there  is  a  need  to  collect  the  existing  information  on 
electric  arcs  for  defining  an  operation  range  in  order  to  in¬ 
vestigate  the  aerodynamic  aspect  of  the  flow  field  phenomenon 
caused  by  electric  discharge. 

The  breakdown  voltage  as  indicated  by  Cobine  (Ref.  18)  and 
explained  in  Appendix  IV  could  be  as  high  as  23,000  volts  for 
the  present  problem.  Of  course,  this  is  not  practical  for  any 
possible  aerospace  applications.  Experimental  methods  for 
establishing  an  arc  with  reasonable  breakdown  voltage  are 
available  for  different  electrode  configurations.  A  literature 
survey  of  the  available  information  under  various  conditions  is 
discussed  in  detail  in  Appendix  IV  and  revealed  the  possibility 
of  establishing  an  arc  at  a  breakdown  voltage  at  about  400 
volts.  The  value  is  encouraging  but  the  information  curves 
are  not  directly  comparable  with  the  present  problem.  Never¬ 
theless,  the  range  of  investigation  can  be  considered  as  a  guide 
for  the  purpose  of  preparing  experimental  investigations. 

The  operating  characteristics  for  a  cylindrical  arc  were 
analyzed  by  Maecker  (Ref.  19)  with  the  assistance  of  some 
experimental  results.  Following  Maecker’ s  approach,  the 
voltage-current  variations  of  an  arc  operating  characteristic 
at  a  nitrogen  stream  of  5000°R  was  evaluated  in  Appendix  V. 

There  the  consideration  of  thermal-chemical  effect  concerning 
ionization  for  an  arc  discharge  at  various  temperatures  and 
pressures  was  based  on  the  results  of  Cambel  (Ref.  20). 
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SECTION  III 


EXPERIMENTAL  PROGRAM 


Electric  arc  studies  have  relied  on  experimental  programs 
because  of  the  complex  emission  processes  caused  by  various 
parameters,  such  as  the  electrode  surfaces  and  the  surrounding 
conditions.  Available  theories  can  only  provide  qualitative 
guidance  for  designing  the  experimental  apparatus.  The  desired 
configuration  with  the  electrodes  embedded  in  the  same  surface 
without  facing  each  other  has  not  been  examined  in  the  litera¬ 
ture.  Preliminary  study  of  Appendix  IV  indicated  that 

i)  the  breakdown  voltage  could  be  as  high  as  23,000  volts 
for  a  pair  of  electrodes  separated  by  one  quarter  of 
an  inch  between  the  shortest  distance,  and 

ii)  the  arc  might  not  be  stable  without  the  confinement 
of  a  strong  magnetic  field. 

Consideration  for  apparatus  design  further  revealed  that 
the  material  problem  could  be  very  serious  especially  in  the 
vicinity  of  the  electric  arc.  In  order  to  study  the  aerodynamic 
phenomena  on  pressure  rises  caused  by  electric  discharge,  all 
the  re1ated  problems  need  to  be  investigated.  Three  different 
types  of  testing  nozzles  were  designed  and  tested.  The 
difficulties  encountered  during  the  testing  processes  are  pre¬ 
sented  in  Appendix  VI . 

3.1  TEST  FACILITY 

The  experimental  program  was  conducted  at  the  Arnold 
Engineering  Development  Center,  Arnold  Air  Force  Station, 
Tennessee.  The  two-megawatt  arc  heater  located  in  the 
Propulsion  Wind  Tunnel  was  selected  because  it  provides  i) 
the  high  temperature  gas  stream  for  simulating  the  rocket  ex¬ 
haust  and  ii)  the  required  d-c  power  supply  for  starting  and 
maintaining  the  electric  arc.  Since  the  facility  has  been 
scheduled  for  other  higher  priority  tests  at  the  same  time, 
the  apparatus  and  instrumentation  were  so  designed  that  a 
minimum  alteration  of  the  standard  set-up  was  needed. 

The  arc  heater  supplied  heated  nitrogen  at  a  stagnation 
temperature  range  between  3000°R  and  6000°R,  and  a  stagnation 
pressure  of  about  50  psia.  With  a  vacuum  system  to  provide  an 
ambient  pressure  approximately  2  psia,  at  a  continuous  flow 
rate  of  .3  Ib/sec,  the  high  temperature  environment  of  a  rocket 
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nozzle  -was  able  to  be  simulated  at  an  exhaust  Mach  number 
approximately  3.  Due  to  the  high  temperature  environment  of  a 
cooling  system  at  400  psia  with  a  maximum  water  flow  rate 
approximately  300  gallons  per  minute  was  used  to  maintain  the 
test  apparatus  within  the  allowable  stress  limit.  A  schematic 
diagram  of  the  test  facility  is  shown  i:i  Fig.  8  and  a  photo¬ 
graph  of  this  facility  is  shown  in  Fig.  9. 

In  order  to  provide  the  electric  energy  for  establishing 
ana  maintaining  an  electric  arc,  a  d-c  power  supply  up  to  1200 
volts  was  used  at  a  maximum  possible  current  of  800  amperes. 

The  d-c  voltage  can  only  be  applied  at  an  increment  of  75  volt 
interval;  the  voltage  and  current  were  registered  on  standard 
meters  and  recorded  photographically  at  three  frames  per  second. 

The  stagnation  temperature  of  the  arc  heated  nitrogen  was 
determined  by  using  the  relation  obtained  from  thermal  energy 
balance  between  the  cooling  water  and  the  heated  nitrogen.  No 
temperature  measurement  was  attempted  in  the  test  section  be¬ 
cause  the  arc  temperature  generally  is  out  of  the  range  of 
conventional  temperature  measurement  technique  such  as  thermo¬ 
couples.  Optical  systems  for  temperature  measurements  were 
studied,  but  the  reliability  and  accuracy  could  not  be  proven 
to  justify  the  effort. 

The  stagnation  pressure  of  the  heated  nitrogen  was  measured 
directly  from  the  stilling  chamber.  The  static  pressures  in 
the  gas  stream  was  measured  through  wall  pressure  taps  by 
transducers  and  manometers,  and  recorded  by  an  oscillograph  and 
cameras  respectively.  Due  to  the  distance  between  the  test 
apparatus  and  the  instrumentation  panel,  a  minimum  time  of  0.8 
seconds  for  the  transducer  readings  and  3.0  seconds  for  the 
manometer  readings  were  required  to  give  the  steady  state 
response.  The  detailed  study  of  the  response  time  is  discussed 
in  Appendix  VI.  The  standard  set-up  of  the  test  facility  is 
described  in  the  AEDC  Test  Facilities  Handbook  (Ref.  22). 

3.2  TEST  NOZZLES 

Three  test  nozzles  were  designed  for  adapting  to  the  two- 
megawatt  heater  for  rocket  nozzle  simulation.  Due  to  the  high 
temperature  requirement,  only  water-cooled  copper  and  ablative 
composite  material  were  chosen  for  this  study. 

3.2.1  Axisymmetric  Copper  Nozzle 

A  water-cooled  axisymmetric  copper  nozzle  was  built  for  the 
purpose  of  studying:  i)  arc  establishment,  ii)  insulation  mate¬ 
rials,  and  iii)  electrode  configuration  investigations.  Prelimi¬ 
nary  information  on  arc  establishment  was  obtained  by  surveying 


the  available  literature  as  outlined  in  Appendix  IV.  The  heat 
transfer  analysis  was  made  by  considering  the  conductive  and 
convective  heat  transfer  at  the  test  conditions.  The  insulation 
materials  were  determined  by  trial  and  error  method  as  outlined 
in  Appendix  VI.  Two  pairs  of  electrodes  were  installed  for 
obtaining  the  information  on  breakdown  voltages  affected  by 
electrode  arrangements.  However,  the  available  space  for 
pressure  measurements  was  limited  by  the  number  of  electrode 
configuration  possibilities.  A  schematic  diagram  of  the 
axisymmetric  copper  nozzle  is  shown  in  Fig.  10;  the  photographs 
of  the  side  and  end  views  are  shown  in  Figs.  11  and  12, 
respectively,  to  indicate  the  pressure  tap  and  electrode 
arrangements.  The  water-cooled  axisymmetric  nozzle  was  modi¬ 
fied  for  studying  pressure  distributions.  A  schematic  diagram 
of  the  modified  axisymmetric  copper  nozzle  is  shown  in  Fig.  13 
to  illustrate  the  pressure  lap  locations  with  respect  to  the 
tested  electrode  configuration.  The  photographs  of  Figs.  14 
and  15  show  the  assembly  and  the  exit  of  this  nozzle. 

3.2.2  Rectangular  Copper  Nozzle 

The  rectangular  copper  nozzle  was  designed  for  the 
following  purposes: 

i)  Visualization  of  electric  arc  and  flow  field  interaction, 

ii)  Electrode  configuration  variations,  and 

iii)  Material  selections  for  electrodes,  insulators,  and 
plastic  nozzle  parts. 

In  order  to  meet  these  requirements,  three  nozzle  sections 
were  used,  namely,  the  throat  section,  the  top  and 
side  wall  sections  and  the  bottom  wall  section.  Water-cooled 
copper  was  used  for  the  first  two  sections  with  quartz  windows 
installed  at  the  downstream  portion  of  the  side  walls  for  flow 
visualization.  Various  materials  were  used  for  the  third 
section  which  was  simply  a  flat  plate  with  electrodes  and 
pressure  taps.  This  third  section  can  be  easily  replaced  for 
investigating  electrode  configurations.  A  drawing  of  the 
rectangular  nozzle  is  shown  in  Fig.  16,  photographs  of  Figs. 

17  and  18  show  the  three  main  sections  of  the  rectangular 
nozzle;  Figs.  19  and  20  show  the  assembly  and  the  side  views, 
respectively . 
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Axisymmetric  Silicon  Phenolic  Nozzle 


The  •water-cooled  copper  nozzle  was  inherently  accompanied 
with  a  cold  boundary  layer  which  may  affect  the  electric  dis¬ 
charge  characteristics.  In  order  to  investigate  the  possible 
effect  on  the  flow  field,  an  axisymmetric  silicon  phenolic 
nozzle  was  built  identical  with  the  axisymmetric  water-cooled 
copper  nozzle.  Due  to  the  high  ablation  rate  of  this  material, 
a  graphite  throat  insert  was  used  for  better  duration.  A 
design  drawing  of  the  axisymmetric  silicon  phenolic  nozzle  is 
shown  in  Fig.  21.  Photographs  of  Figs.  22  and  23  show  the 
pressure  tap  locations  and  overall  assembly  of  the  nozzle. 

3.3  TEST  PROCEDURE 

After  all  instruments  were  properly  calibrated,  the  test 
was  started  for  any  predetermined  test  conditions  .  The  standard 
procedure  for  each  test  run  was  as  follows: 

The  vacuum  system  was  started  first  to  establish  the 
desired  back  pressure  at  the  nozzle  exit .  The  cooling  water 
was  then,  started  to  circulate.  The  water  flow  rates  of  all 
cooling  sections  were  recorded  and  inspected  to  insure  sufficient 
amount  of  cooling  for  preventing  any  possible  overheating  of 
any  parts.  Since  the  test  was  conducted  at  very  high  temperature 
and  low  pressure  environment,  the  cooling  water  flow  rate  could 
be  very  critical  in  order  to  keep  the  apparatus  within  the 
allowable  stress  limit. 

When  the  cooling  water  system  was  functioning  properly, 
the  cold  nitrogen  flow  was  then  started.  Simultaneously,  all 
the  cameras  and  oscillograph  were  turned  on  to  record  time 
dependence  variations  of  all  pressures  and  temperatures . 

The  arc  heater  was  then  started  to  heat  the  nitrogen  at 
the  stilling  chamber.  As  soon  as  steady  state  was  reached 
(generally,  it  required  3  to  5  seconds),  the  d-c  power  for 
electric  discharge  was  applied,  A  high  frequency  starter  was 
connected  in  the  circuit  and  cculd  be  used  if  necessary. 

The  operating  characteristics  could  be  established  by 
decreasing  the  circuit  resistance  or  increasing  the  circuit 
voltage.  The  variations  of  both  current  and  voltage  were 
continuously  recorded  by  photographing  the  ammeter  and  volt¬ 
meter  readings.  The  voltage  for  electric  discharge  was  dis¬ 
continued  after  a  period  of  time  which  was  considered  long 
enough  for  steady  response  of  all  the  pressure  measurements. 

The  nitrogen  flow,  instrumentation  recorders,  cooling  water 


and  the  vacuum  system  were  then  turned  off  in  sequence.  In 
some  cases,  the  static  pressure  after  discharge  did  not  return 
to  the  before-discharge  readings  as  discussed  in  Appendix  VI . 
Cautions  were  taken  to  allow  at  least  3  seconds  between  the 
voltage  off  and  instrumentation  off  in  order  to  check  the 
possible  cause  of  pressure  rise  other  than  the  electric  dis¬ 
charge  .  The  actual  pressure  rise  caused  by  electric  discharge 
was  determined  by  using  the  pressure  difference  between  the 
pressure  reading  during  and  after  discharge.  This  method  of 
determining  pressure  rise  caused  by  electric  arc  is  being 
considered  as  conservative  because  the  material  is  deforming 
continuously  even  shortly  after  the  electric  discharge. 

3.4  RESULTS  AND  DISCUSSION 

3.4.1  Breakdown  Voltage  in  a  Supersonic  Stream 


No  information  is  available  concerning  the  breakdown  voltage 
between  two  parallel  electrodes  embedded  flush  in  the  same 
surface  in  a  supersonic  stream.  Various  sources  as  discussed 
in  Appendix  IV  indicated  the  possible  breakdown  voltage  ranges 
from  400  to  23,000  volts  for  the  desired  configurations.  This 
uncertainty  caused  a  great  deal  of  difficulty  in  selecting  an 
appropriate  test  facility  for  the  experimental  investigation. 

A  compromise  had  to  be  made  in  order  to  cover  the  reasonable 
ranges  for  the  flow  Mach  numbers,  the  stream  temperatures,  and 
the  d-c  power  availability.  The  two-megawatt  heater  provided 
the  conditions  closest  to  the  desired  range  with  the  disad¬ 
vantages  that  the  flow  Mach  number  was  limited  to  3  and  the 
side  force  could  not  be  directly  measured. 

The  axisymmetrical  water-cooled  copper  nozzle  was  used  for 
investigating  the  breakdown  voltage.  A  75  volt  increment  was 
applied  for  testing  at  various  design  parameters.  The 
experimental  results  are  tabulated  in  Table  n.  It  can  be 
seen  that  a  breakdown  voltage  of  375  volts  was  needed  to  start 
an  arc  between  a  pair  of  electrodes  of  one-quarter  inch  apart 
mounted  flush  on  the  same  surface  over  which  a  Mach  2.5 
nitrogen  stream  flowing  at  a  static  pressure  of  2  psia  and  a 
stagnation  temperature  of  4000°R.  The  breakdown  voltage  was 
substantially  reduced  by  increasing  the  stream  temperature, 
decreasing  the  distance  between  the  electrodes,  arranging  the 
electrodes  parallel  instead  of  perpendicular  to  the  flow  direction, 
or  seeding  the  stream  with  low  ionization  potential  material. 

It  is  to  be  emphasized  that  the  tabulated  results  in  Table  II 
are  only  qualitative  nature  because  of  the  d-c  power  supply  was 
limited  to  75  volt  increments.  Nevertheless,  the  possibility 
of  successfully  establishing  an  electric  arc  in  a  supersonic 
stream  is  proved  to  be  within  the  practical  voltage  range > 
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TABLE  II  -  ARC  ESTABLISHMENT  PARAMETERS 
(MACH  NO.  OF  GAS  STREAM  =2.5) 


Investigated  Parameters 

Circuit 

Voltage 

(Vo) 

Total 

Temp. 

(Tt) 

Higi  Freq. 
Starter 
(S) 

Dist. (d) 
Between 
Electrode 

Electrode 
Arrg ' t 

Polarity 
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(K2Co3) 

°R 

— 

in. 

— 

— 

% 
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4200 
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.25 
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300 

4200 
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.125 
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ll 

— 

0 

150 

Reverse 
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w 

.20 

ll 

Reverse 

0 

225 

.25 
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1  1 

5700 


w 


20 


28 


75 


namely,  the  electric  field  required  to  initiate  discharge  was 
about  2  x  103  volts  per  inch  at  the  studied  condition.  More¬ 
over,  the  use  of  high  frequency  starter  was  found  to  be  very 
effective  in  assistance  to  establish  the  arc. 

3.4.2  Operating  Characteristics  of  an  Electric  Arc 

The  operating  characteristics  of  an  electric  arc  can  only 
be  analyzed  for  highly  simplified  models  as  discussed  in  Appendix 
V.  Conf igurations  for  the  proposed  scheme  have  not  yet  been 
studied.  The  experiment  for  establishing  a  voltage-current 
curve  was  conducted  by 

i)  increasing  the  circuit  voJtage,  V  ,  as  shown  in 
Figure  24  and,  ° 

ii)  reducing  the  circuit  resistance,  R  . 

c 

The  operating  characteristics  are  shown  in  Figure  25  for  the 
electrode  arrangements  perpendicular  and  parallel  to  the  flow 
direction. 

3.4.3  Pressure  Distributions 


The  pressure  distributions  in  the  vicinity  of  the  elec¬ 
trodes  were  measured.  However,  due  to  the  high  temperature 
environment  (approximately  10,000°K  arc  temperature),  and  the 
high  velocity  (about  5000  ft/sec)  flow  field,  the  material 
problem  became  so  severe  that  the  pressure  distributions 
could  not  be  obtained  without  any  interference  from  the  material 
deformation.  Three  test  nozzle  designs  were  used  for  in¬ 
vestigating  the  static  pressure  distributions: 

3. 4. 3.1  Axisyrametric  Copper  Nozzle 

The  pressure  distribution  along  the  nozzle  centerline  is 
shown  in  Figs.  26,  27,  and  28  for  the  electrode  configurations, 

A,  B,  and  C,  respectively.  Configuration  A  was  for  electrodes 
one-quarter  inch  apart  perpendicular  to  the  flow  direction. 
Configuration  B  was  for  electrodes  one-eighth  inch  apart  perpen¬ 
dicular  to  the  flow  direction.  Configuration  C  was  for 
electrode  one-quarter  inch  apart  parallel  to  the  flow  direction. 

The  area  affected  by  pressure  rise  in  configuration  A  is  compara¬ 
tively  larger  than  those  of  configurations  B  and  C.  However, 
the  pressure  taps  were  not  dense  enough  to  provide  the  detailed 
distribution  especially  in  the  vicinity  of  the  electrodes.  Further¬ 
more,  the  deformation  of  the  insulation  material  destroyed  the 
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smooth  contour  of  the  nozzle  surface,  which  caused  the  static 
pressures  after  discharge  to  differ  from  that  before  discharge. 
The  insulation  damage  can  be  seen  from  a  photograph  for  electrode 
after  discharge  in  Fig.  2d.  By  examining  the  post  discharged 
test  nozzle  the  different  colors  on  the  wall  downstream  of  the 
electrodes  indicated  a  very  similar  pattern  as  observed  down¬ 
stream  of  the  injection  part  of  the  secondary  fluid  injection 
cases  (Ref.  23).  These  color  changes  were  considered  as  an 
evidence  of  the  presence  of  a  shock  wave  caused  by  the  mass 
addition  in  the  secondary  fluid  injection  case  and  by  the 
electric  arc  discharge  in  the  present  case . 

The  insulation  material  created  a  serious  problem,  as 
discussed  in  Appendix  VI.  Among  all  the  tested  materials,  boron 
nitride  appeared  to  be  the  only  one  possible  to  provide  the 
desired  performance  and  was  used  for  electrode  insulators  for 
all  the  remaining  experiments.  The  electrode  configuration  C 
was  tested  for  measuring  pressure  distributions  at  various  power 
inputs.  The  arc  was  established  at  a  breakdown  voltage  of  75 
volts,  with  the  assistance  of  a  high  frequency  starter  under  the 
flow  conditions  of  a  Mach  2.5  nitrogen  stream  seeded  with 
K2C03  (approximately  0.28%  of  potassium)  at  5700°R  stagnation 
temperature  and  2  psia  static  pressure.  The  pressure  distribu¬ 
tions  are  shown  in  Fig.  30.  The  history  of  the  pressure 
variations  corresponding  to  the  power  variations  is  shown  in 
Fig.  31.  It  can  be  seen  that  the  pressure  rise  is  a  function 
of  power  input.  Moreover,  the  static  pressure  after  discharge 
returned  to  the  values  of  before  discharge  which  proved  that 
the  boron  nitride  could  be  used  as  insulators  between  the 
electrodes  and  the  copper  wall. 

It  is  necessary  to  point  out  that  the  static  pressure  at 
the  nozzle  exit  was  not  following  the  isentropic  expansion  curve 
even  without  electric  discharge.  This  was  caused  by  flow 
separations  resulting  from  the  leaking  problem  of  the  vacuum 
system  which  prevented  the  test  facility  from  reaching  the 
designed  nozzle  back  pressure.  The  axisymmetric  copper  nozzle 
was  modified  by  changing  the  exit  area  to  the  attainable  back 
pressure  of  the  facility  and  installing  more  pressure  taps  in 
the  vicinity  of  the  electrodes.  The  problems  of  flow  separations 
and  insulation  damages  appeared  to  be  considerably  improved. 

The  static  pressure  distributions  are  shown  in  Fig.  32.  The 
static  pressure  rise  due  to  electric  arc  discharge  seemed  to 
affect  only  a  very  small  area. 
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3. 4. 3. 2  Rectang'i  3  ar  Copper  Nozzle 

The  pressure  distribution  of  the  rectangular  copper 
nozzle  was  measu,'(,d  using  several  types  of  material  as  the 
bottom  nozzle  wa  .  The  purpose  of  eliminating  the  cold 
boundary  layer  e:i  j  ect  was  the  main  reason  for  using  a  non- 
metallic  material  as  the  bottom  nozzle  .wall.  However,  several 
problems  were  encountered  during  the  testing.  Boron  nitride 
was  the  first  material  selected  because  of  its  high  thermal 
shock  resistance,  but  the  difficulty  of  preventing  leaks 
between  the  pressure  taps  and  the  boron  nitride  plate  made 
this  material  impractical  from  the  cost  standpoint.  The  leak 
prevention  effort  for  using  boron  nitride  is  discussed  briefly 
in  Appendix  VI . 

Silicon  phenolic  material  was  used  with  fair  success 
when  the  electrodes  were  insulated  with  boron  nitride.  The 
insulation  as  necessary  because  charred  silicon  phenolic  is  a 
good  electric  conductor.  The  pressure  distributions  are  shown 
in  Fig.  33  for  electrode  arrangement  perpendicular  to  the  flow 
direction.  The  flow  field  disturbances  caused  by  electric 
arc  were  recorded  with  high  speed  movie  cameras  for  visualization 
purposes  (available  on  request).  Fig.  34  shows  a  typical 
movie  frame  of  the  flow  field  disturbance  due  to  electric 
discharge  which  is  very  similar  to  the  picture  for  secondary 
fluid  injection  as  shown  in  Ref.  3.  It  was  observed  that 
the  extensive  heat  of  the  arc  caused  the  plate  to  ablate  and 
deform  and  prevented  the  pressure  reading  from  returning  to 
the  original  value  after  discharge.  Moreover,  the  highest  p^ak 
of  the  pressure  rise  was  not  detected  because  the  insulation 
occupied  the  space  where  the  pressure  tap  was  located  in  the 
axisymmetric  nozzle.  The  phenolic  plate  after  discharge  is  shown 
in  Fig.  35.  A  water-cooled  copper  plate  was  also  used  for  the 
configuration  of  electrodes  perpendicular  to  the  flow  direction. 
The  electric  discharge  was  established  through  the  copper  plate. 
Ceramic  coating  on  the  surface  of  the  water-cooled  copper  plate 
was  used.  The  pressure  distribution  is  shown  in  Fig.  36.  The 
static  pressure  returned  to  the  original  values  after  discharge. 
The  ceramic  coated  copper  plate  seemed  to  be  a  better  approach 
to  prevent  erosion  as  well  as  to  maintain  the  arc  between 
electrodes.  However,  the  ceramic  coating  could  not  remain 
effective  for  sufficient  time  to  give  steady  state  manometer 
readings  and  the  effective  pressure  rise  area  appeared  to  be 
very  small.  The  picture  in  Fig.  37  shows  the  copper  plate 
after  discharge. 

It  is  noted  that  the  supplied  electric  energy  partially 
went  to  ablate  the  material  as  in  the  silicon  phenolic  plate 
case  and  to  heat  the  cooling  water  as  the  ceramic  coated  copper 
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plate.  In  any  event,  the  phenomenon  that  electric  discharge 
can  cause  pressure  rise  in  a  supersonic  flow  field  is  confirmed. 

The  effective  area,  however,  is  very  snail  for  free  stream 
Mach  numbers  near  3. 

Ac  electrode  configuration  arranged  parallel  to  the  flow 
direction  was  also  tested.  No  substantial  pressure  rise  was 
measured  during  discharge  for  all  the  pressure  taps.  The  ceramic 
coated  copper  plate  after  discharge  is  shown  in  Fig.  38. 

3. 4. 3. 3  Silicon  Phenolic  Axisymmetric  Nozzle 

The  pressure  distributions  obtained  from  the  silicon 
phenolic  axisymmetric  nozzle  are  shown  in  Fig.  39.  The  static 
pressure  after  discharge  failed  to  return  to  the  before  discharge 
values  because  of  the  severe  damage  of  the  surface  of  phenolic 
material.  The  information  on  pressure  rise  was  obtained  by 
using  the  pressure  differences  between,  during  and  after  dis¬ 
charges.  The  actual  pressure  rise  caused  by  electric  discharge 
could  be  higher  than  the  used  value  because  of  the  continuous 
deformation  of  the  phenolic  material.  Nevertheless,  this 
method  provided  the  quantitative  results  on  the  most  conservative 
side.  Thus,  the  concept  that  the  electric  discharge  can  cause 
sudden  pressure  rise  is  confirmed.  Fig.  40  shows  the  picture 
of  the  phenolic  nozzle  taken  after  discharge.  By  examining 
the  post  discharged  test  nozzle,  the  color  changes  on  the  wall 
downstream  of  the  electrodes  indicates  the  shock  wave  pattern 
caused  by  the  electric  arc. 

3.5  COMPARISON  OF  PRESSURE  RISE  IN  MAIN  STREAM  DUE  TO  ELECTRIC 

DISCHARGE  AND  THE  PRESENCE  OF  SOLID  BODY 

It  has  been  reported  by  Roman  (Ref.  17)  that  a  magnetic  confined 
arc  in  a  cross  flow  behaves  like  a  heated  solid  body  as  has  been 
mentioned  in  Section  II.  In  the  present  experiment,  the  arc 
discharge  is  not  magnetically  confined.  However,  it  has  been 
observed  that  the  same  arc  -  solid  body  analogy  still  exists. 

Once  the  discharge  is  established,  energy  is  continuously  being 
fed  in.  The  discharge  takes  a  curved  path  rather  than  a  straight 
one  (as  in  the  magnetically  confined  arc)  because  of  the  blowing 
effect  of  the  main  stream. 

Figure  41  shows  the  pressure  rise  caused  by  the  protrusion 
of  electrodes  above  the  nozzle  wall  about  0.06  in.  without  discharge 
and  caused  by  electric  discharge  with  flush  mounted  electrodes. 

It  is  clear  that  the  pressure  rise  caused  by  these  two  effects  are 
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quite  similar.  The  pressure  rise  was  shifted  downstream  in  the 
discharge  cases  (Fig.  41).  This  was  caused  by  the  curved  shape 
of  the  discharge  due  to  the  blowing  effect.  Vas  and  Bogdonoff 
(Ref.  24)  reported  that  the  pressure  just  upstream  of  a  forward 
facing  step  increases  linearly  with  step  height  as  long  as  the 
step  height  is  less  than  the  boundary  layer  thickness.  Since 
the  distance  of  protrusion  of  electrodes  were  less  than  the 
boundary  layer  thickness  (about  0.2  in.  in  our  test  case), 
it  may  be  able  for  us  to  estimate  the  equivalent  step  height 
corresponding  to  a  similar  pressure  rise  by  an  arc  discharge. 

3.6  SIDE  FORCE  EVALUATION 

The  side  force  was  evaluated  by  integrating  the  pressure 
differences  during  and  after  discharge  over  the  affected  areas. 
The  detailed  evaluation  procedures  are  discussed  in  Appendix  VII. 
Due  to  the  limited  experimental  data  available,  this  evaluation 
only  gives  the  order  of  magnitude  of  the  generated  side  force. 

The  evaluated  results  are  compared  with  theoretical  calculations 
and  presented  in  Table  III.  The  experimental  level  for  the 
parameter  Fs/Pin  is  low  but  comparable  to  that  of  the  theoretical 
calculations.  The  energy  lost  in  ablating  electrodes,  insulators, 
and  conduction  to  the  cooled  walls  has  not  been  considered. 


31 


TABLE  III.  COMPARISON  OF  MEASURED  AND  PREDICTED  SIDE  FORCE 
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SECTION  IV 


CONCLUSIONS 


Analytical  and  experimental  investigations  were  conducted 
for  the  electric  discharge  in  a  supersonic  stream.  The 
following  conclusions  were  reached. 

4.1  STATIC  PRESSURE  RISE  CAN  BE  PRODUCED  BY  ELECTRIC  DISCHARGE 

A  static  pressure  rise  occurred  in  the  vicinity  down¬ 
stream  of  the  electrodes  when  an  arc  discharge  was  applied  in 
a  supersonic  flow  field.  This  type  of  static  pressure  distribu¬ 
tion  was  also  observed  in  secondary  fluid  injection  case, 
especially  when  gas  was  used  as  the  inject ant.  The  pressure 
peak  caused  by  arc  discharge  agreed  with  the  trend  predicted 
by  the  blast  wave  theory. 

4.2  SIDE  FORCE  GENERATED  BY  AN  ARC  DISCHARGE  MAY  BE  PREDICTED 

WITH  SIMPLIFIED  THEORY 

The  blast  wave  analogy,  linearized  small  disturbance  theory, 
and  plane  flame  theory  all  predictec  with  the  same  order  of 
magnitude  of  the  experimental  result.  The  blast  wave  analogy 
appeared  to  be  more  advantageous  because  it  considered  the 
nonlinear  interaction  phenomena  of  the  energy  source  with  the 
flow  and  also  resulted  in  a  simple  mathematical  expression. 

4.3  ARC  CAN  BE  ESTABLISHED  AND  MAINTAINED  IN  THE  SUPERSONIC 

STREAM 

There  is  no  report  available  on  the  arc  discnarge  to  a 
pair  of  electrodes  which  are  embedded  flush  on  the  wall  and 
are  not  directly  facing  each  other.  From  this  study,  it  was 
found  that  the  arc  can  be  established  and  maintained  in  a 
supersonic  flow  without  magnetic  confinement. 

4.4  THE  BREAKDOWN  VOLTAGE  IS  WITHIN  THE  PRACTICAL  RANGE 

The  electric  field,  to  initiate  arc  discharges,  was 
about  2  x  103  volts/in.  for  nitrogen  stream  at  Mach  number  2.5 
with  static  pressure  of  2  psia  and  stagnation  temperature  of 
4000°R.  The  required  electric  field  can  be  reduced  to  approxi¬ 
mately  3  x  102  volts/in.,  if  the  gas  stream  stagnation 
temperature  is  about  6000°R  with  0.28%  of  potassium  seeds  and 
favorable  electrodes  orientation. 
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4.5  THE  SIDE  FORCE  SPECIFIC  IMPULSE  OF  ARC  DISCHARGE  DEPENDS 

ON  THE  SPECIFIC  ENERGY  OF  POWER  SOURCE 

Comparing  the  specific  impulse  of  arc  discharge  scheme 
to  that  of  the  secondary  fluid  injection,  method,  it  was 
found  that  the  arc  discharge  technique  only  becomes  competitive 
when 

i)  high  specific  energy  power  source  is  available 
for  rocket  type  vehicles,  or 

ii)  sufficient  power  source  is  available  on  board  for 
re-entry  type  vehicles. 

4.6  MATERIAL  PROBLEMS  ASSOCIATED  WITH  ELECTRIC  ARC  DISCHARGE 

IN  SUPERSONIC  FLOW 

It  was  found  in  the  experiments  that  the  available 
material  can  only  sustain  a  very  short  duration  without  serious 
deformation.  In  order  to  obtain  accurate  measurements,  highly 
sensitive  instrumentation  must  be  used.  The  experimental 
program  was  discontinued  because  the  development  of  high 
temperature  material  at  the  present  is  not  adequate  for 
system  applications. 


SECTION  V 


RECOMMENDATIONS 


Based  on  the  results  of  this  investigation,  the  following 
recommendations  can  be  made: 

5.1  ARC  DISCHARGE  FOR  ROCKET  APPLICATIONS 

The  most  attractive  place  to  use  the  arc  discharge  tech¬ 
nique  is  for  nuclear  energy  propelled  rockets  where  energy 
consumption  for  auxiliary  power  does  not  result  in  significant 
weight  and  cost  penalties.  For  chemical  rockets,  however,  the 
only  possibility  of  making  the  electric  arc  discharge  for 
thrust  vector  control  to  be  feasible  is  to  employ  high  specific 
energy  power  source  which  may  be  available  in  the  future. 

5.2  ARC  DISCHARGE  FOR  RE-ENTRY  APPLICATIONS 

For  the  long  time  mission  space  vehicles,  the  solar  energy 
may  be  stored  and  used  for  the  re-entry  trajectory  control 
purpose.  The  phenomena  of  arc  discharge  in  a  very  high  Mach 
number  case  needs  further  investigation. 

5.3  A  COMBINED  APPLICATION  OF  SECONDARY  MSS  INJECTION  AND 

ARC  DISCHARGE 

From  the  side  force  specific  impulse  point  of  view,  it  is 
found  that  the  arc  discharge  alone,  or  purely  by  the  thermal 
energy  addition,  under  the  present  available  power  source  is 
not  as  effective  as  the  secondary  injection.  A  new  device  may 
be  developed  by  combining  the  two  methods .  The  arc  discharge 
could  be  designed,  such  as  the  arc  jet  device,  to  provide  thermal 
energy  for  the  secondary  injection  gas  rather  than  directly 
delivered  to  the  main  gas  stream.  By  this  means,  the  secondary 
gas  flow  will  receive  an  additional  energy  or  to  increase  its 
momentum  as  it  leaves  the  injection  port,  and  resulted  in  a 
mere  effective  side  force  generation  than  by  either  gas  injection 
or  arc  discharge  alone . 

5.4  ARC  APPLICATIONS  FOR  SUPERSONIC  COMBUSTION 

The  electric  arc  discharge  proved  to  be  a  very  practical 
means  for  energy  addition  at  low  pressure  environment.  This 
phenomenon  may  assist  to  overcome  the  difficulties  for  diffussional 
mixing  in  supersonic  combustion  at  low  pressure.  There,  in 
general,  the  long  relaxation  time  results  in  a  combustor  being 
too  long  for  practical  purposes.  Arc  discharge  in  the  fuel  and 
air  mixing  region  may  accelerate  the  mixing  and  combustion 
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processes.  Moreover,  the  required  power  source  can  be  obtained 
through  a  MHD  generator  which  may  be  located  at  the  downstream 
portion  of  the  system. 
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Figure  1  .  Flight  Direction  Control  By  Electric 

Discharge 


Electrodes 

SIDE-VIEW 


Figure  2.  Schematic  Drawing  of  Generated 
Shock  Wave  Due  to  Arc  Discharge 
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Figure  7.  Theoretical  Performance  of 
Side  Force  Specific  Impulse 
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FIGURE  12.  END-VIEW  OF  THE  AXISYMMETRIC  WATER-COOLED  COPPER  NOZZi.E 


2°  Outlet 


Fig-  14.  Assembly  of  the  Modified  Axisymmetric  Water-Cooled  Copper  Nozzle 


Figure  16.  DRAWING  OF  RECTANGULAR  TEST  NOZZLE 
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Fig.  19.  End-View  of  the  Rectangular  Nozzle  Assembly  with  Silicon 
Phenolic  Plato 


Fig.  23.  End-View  of  the  Axisymmetric  Silicon  Phenolic 


SYMBOLS 


O  =  Axisymmetric  water-cooled  copper  nozzle 
A  =  Rectangular  nozzle  with  copper  plate 
□  =  Modified  axisymmetric  water-cooled  copper  nozzle 
B  =  Axisymmetric  silicon  phenolic  nozzle 
&  =  Rectangular  nozzle  with  silicon  phenolic  plate 
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Figure  25.  Measured  Arc  Characteristic 
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Figure  32 . 


(Rectangular  Nozzle  with  Silicon  Phenolic  Plate 
Cross-Flow,  0.25  in. apart)  ‘ 
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Electric  Arc  ana  Flow  Field  Interaction 


Figure  34.  Flow  Field  Disturbance  Visualization 


FIGURE  35  POST-DISCHARGE  OF  THE  RECTANGULAR  NOZZLE (CROSS-FLOW 

SILICON  PHENOLIC  PLATE) 


FIGURE  37  POST-DISCHARGE  OF  THE  RECTANGULAR  NOZZLE (CROSS-FLOW 

COPPER  PLATE) 


FIGURE  38  POST-DISCHARGE  OF  THE  RECTANGULAR  NOZZLE (PARALLEL- 

COPPER  PLATE) 


FIGURE  40.  POST-DISCHARGE  OF  THE  AXISYMMETRIC  SILICON  PHENOLIC  NOZZLE 
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APPENDIX  I 


DISCUSSION  ON  "ANALYSIS  OF  THE  FLUID  MECHANICS  OF 
SECONDARY  INJECTION  FOR  THRUST  VECTOR  CONTROL," 

BY  BROADWELL 


Broadwell ~4Ref  * — i)  demonstrated  that  the  blast  wave  analogy 
can  be  employed  in  explaining  and  predicting  the  side  force  by 
the  secondary  mass  injection,  thrust  vector  control.  He  suc¬ 
cessfully  derived  the  side  force  expression  by  integrating 
the  pressure  differential,  obtained  by  the  blast  wave  theory, 
over  the  affected  area  and  resulted  in 


F„  =  cr( 7)  E  M 

s  *  00 
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Obviously,  to  our  study,  the  line  energy  source  E  is 
contributed  by  the  rate  of  thermal  energy  addition  Q  or  equiv¬ 
alent  to  the  electrical  power  input  PiQ  only.  The  relation 
of  E  to  P-n  has  been  discussed  in  detail  in  the  main  text. 

In  case  oisecondary  mass  injection,  Broadwell  related  E  to 
the  drag  force  caused  by  the  secondary  fluid  injection  which 
is  analogous  to  the  derivation  originated  in  Lees'  paper 
(Ref .  16) .  Lees  argued  that  in  order  to  hold  a  blunt-nosed 
body  in  a  hypersonic  stream,  it  is  necessary  to  apply  a  drag 
force.  This  drag  force  then  can  be  shown  equivalent  to  the 
line  energy  source  in  the  cylindrical  blast  wave  case. 

The  linearized  small  disturbance  theory  was  employed  by 
Wu,  Chapkis  and  Mager  (Ref.  2)  in  determining  the  side  force 
generated  by  the  slot  type  secondary  fluid  injection.  The 
theoretical  prediction  yields  good  agreement  with  experimental 
results.  However,  Broadwell  (Ref.  4)  stated  that  the  linearized 
disturbance  theory  implicitly  included  the  mass  source  entering 
the  free  stream  with  an  x-component  velocity  of  the  same 
magnitude  of  the  free  stream.  Then  he  argued  that  the  initial 
x- velocity  of  inject ant  has  to  be  zero,  and  added  an  additional 
term  to  the  side  force  equation  originally  derived  in  Wu,  et  al, 
(Ref.  2).  His  result  is  as  follows: 
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As  it  nay  be  seen,  and  also  pointed  out  by  himself,  that  the 
added  second  tern  could  be  order  of  magnitude  larger  than 
the  original  first  term,  since  >  1.  Experimental  results 
contradict  with  this  order  of  magnitude  large  side  force 
prediction.  Furthermore,  he  claimed  that  when  »  1  the 
linearized  theory  with  the  added  term  are  comparable  with  the 
blast  wave  theory  with  a  volume  effect.  These  two  theories 
are  not  comparable,  because  the  linearized  theory  is  valid 
only  at  moderate  Mach  number,  while  the  blast  wave  theory  is 
valid  only  at  hypersonic  conditions. 

Ir  fact,  the  mass  enters  the  main  stream  with  zero 
x-component  of  velocity  relative  to  the  surface.  This  mass 
injection  disturbs  the  main  flow.  It  is  a  false  statement 
to  assume  the  flow  enters  with  x-component  velocity  of  u  , 
and  then  ’'make*'  the  x-component  to  be  zero.  The  small 
disturbance  theory  is  derived  under  the  assumption  that  the 
uniform  velocity  is  much  larger  than  the  disturbance  velocity. 
This  is  somewhat  different  from  the  blast  wave  theory, 
because  the  blast  wave  theory  was  originally  derived  from  the 
case  of  a  quiescent  ambient  atmosphere.  Therefore,  in  order 
to  use  the  blast  wave  theory  in  a  steady  flow  stream,  a 
superposition  of  the  same  magnitude  of  the  free  stream  velocity 
to  the  disturbance  source  is  required. 

Therefore,  it  can  be  concluded  that  the  second  term  of 
the  right-hand  side  of  Equation  2  should  not  be  included. 
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APPENDIX  II 


EILLMARTH'S  HEAT  AHD  MASS  ANALOGY 


In  this  Appendix  we  are  going  to  outline  briefly  Willnarth 
analysis  of  the  heat  source,  rsass  source  and  body  force 
analogy  in  supersonic  flow  by  snail  disturbance  theory. 

The  governing  equations  of  notion,  including  nass  source, 
body  force  and  heat  addition,  neglecting  viscosity  and  heat 
conduction,  can  be  written  as  follows: 


7  •  (pW)  =  n  (x,y,z,t)  (1) 

St 

*> 

=  -  Vp  -r  f  (x,y,z,t)  (2) 

at 


DS 

Dt 


CV(Y-D  . 

-  q  (x,y,z,t) 

P 


(3) 


-> 


where 


And,  in 


W  =  velocity  vector 
S  =  entropy  of  the  flow 

h  =  fluid  source  strength,  mass/sec/volume 

"7 

f  =  body  force/volume 

q  =  heat  source  strength,  ft-lb/sec/volume 
addition,  we  assume  the  equation  of  state 


p  =  pRT 


(4) 
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and  the  equation  of  entropy 


S  =  C  In 
v 


(5) 


If  the  disturbances  are  snail,  nanely. 


S  =  S  4-  S 

CO  i 


« 1 

«  1 

« 1 

« 1 

« 1 


Where  W  »  disturbed  velocity  vector  =  uix  4-  viy  4-  wiz.  There 
fore,  the  governing  equations  may  be  linearized.  The  results 
are: 


&p  -* 

u  — i-  4-  pM  V  '  W  =  n  (6) 

dX  x 


82 


(7) 


-> 


3w  -> 


Combining  the  above  equations,  the  result  in  the  x-conponent 
becomes:  *. 


(1  _  H=)  ^  +  21  =  A.  +  2=1.  4  -  J"_  1  (10) 

©  ^  'v  O  O  X 

ox  oy  Sy  p„  p„aa  Ptoa» 


Which  Is  the  fora,  we  cited  in  Section  II  of  Equation  15. 

It  can  be  seen  that  a  mass  source  per  volume  of  strength 

n  will  produce. the  sane  disturbance  as  an  energy  source  of 

strength  q  if  n  =  (7-I)  q/a^,  or  the  sane  as  a  body  force  per 

volume  f  ,  if  n  =  u^f^/a^.  Therefore,  there  exists  an  analogy 

among  the  three  disturbance  sources.  The  equivalences  for 

the  total  rate  of  mass  addition  m,  energy  addition  Q  and  the 

total  body  force  F  are  as  follows: 

x 


00 


(ID 
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APPENDIX  III 


SDR VET  ON  CHEMICAL  BATTERY  POWER  SOURCES 


At  present,  the  available  chemical  battery  power  sources 
are  unsuited  for  use  in  rocket  application.  This  .ls  largely 
due  to  the  fact  that  the  available  batteries  are  not  designed 
for  such  purposes  because  of  its  weight  and  life  duration. 

The  theoretical  performances  of  the  state  of  art  of 
available  chemical  batteries  are  listed  as  follows  (Ref .  25) : 


Type  of  Battery 

Theoretical  Maximum 
Performance,  w-hr/lb 

Lead-acid 

75 

Nickel-cadium 

90 

Mercury-zinc 

100 

Zinc-silver  oxide 

180 

The  actual  performance  of  the  above  batteries  are  considerably 
less  than  the  theoretical  values.  The  extensive  research 
now  under  way  has  substantially  increased  the  battery  performance, 
yet  from  the  system  point  of  view,  it  is  still  too  heavy  for 
the  airborne  type  application. 

Ray  and  Ross  (Ref.  26)  reported  that  a  system  producing 
as  high  as  ^300  w-hr/lb  specific  energy  battery  is  being 
developed.  Its  reactants  are  liquid  lithium  and  gaseous 
chlorine,  and  the  electrolyte  is  fused  lithium  chloride. 
Experimental  data  indicates  a  power  density  of  about  15  w/sq. 
cm.  They  reported  that  the  major  uncertainty  of  this  battery 
is  lifetime.  The  short  lifetime  may  not  be  suited  for  some 
other  purpose  but  could  be  suited  for  the  airborne  short  time 
application. 

It  is  worthwhile  to  point  out  here  that  the  airborne  MHD 
power  source,  besides  the  chemical  battery,  may  be  also  promising 
in  the  future.  J.  B.  Dicks  and  Associates  (Ref.  27)  are  now 
engaged  in  developing  MHD  generators  by  using  superconducting 
magnet  for  airborne  purposes.  Using  a  liquid  propellant 
motor  with  mass  flow  rate  of  4  Ib/sec,  with  the  power  generation 
unit  located  at  the  nozzle  exit  at  a  pressure  of  5  atm,  1*5 
kw/lb  can  be  generated  with  a  maximum  voltage  of  3000  volts . 

This  weight  includes  the  magnet,  channel  and  magnet  coolant. 
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In  order  to  get  some  feeling  of  the  required  power  source 
weight,  a  sample  calculation  for  a  main  thrust  of  P0  =  10,000 
lb  motor  of  chamber  temperature  of  T-fc  =  5000°R,  has  been 
made.  Assuming  that  the  deflection  angle  e  is  4°  and  the 
electrodes  are  embedded  at  Mw  =  4.0  station  inside  the 
exhaust  nozzle,  then  the  required  power  source  weight  consump¬ 
tion  rate  w  can  be  calculated  by  the  following  equation: 


F  tan  e  ^/yR 
[ - FF7^1 - 3 


where  R  =  gas  constant,  o( y)  is  a  constant  from  blast  wave 
theory  and  Esp  ^specific  energy  of  power  source.  Prom  the 
above  equation  the  total  weight  of  power  supply  W  is  there¬ 
fore  directly  proportional  to  the  specific  energy  of  power 
source  E  and  the  duration  of  time  applied.  At  present, 
the  best  battery  available  on  market  (Ref.  28)  is  about  in 
the  order  of  E§p  ~  30  (w-h^/lb.  The ^ corresponding  battery 

consumption  rate  can  be  estimated  as  w  110  lb/sec .  This 
is  still  under  the  assumption  that  a  very  short  time  battery 
is  available.  However,  according  to  the  manufacturer,  the 
battery  is  designed  for  three  minutes  use.  If  the  battery 
of  specific  energy  of  ^  300  (w~hr)/lb  may  be  available  in 
the  future  as  reported  by  Ray  and  Ross  (Ref .  26) ,  the  weight 
consumption  rate  will  be  about  w  ~  11  lb/sec  for  this  specific 
application. 
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iPPENBIX  IV 


INFLUENCE  PARAMETERS  ON  TEE  ARC  DISCEARGE 


In  order  to  determine  the  arc  voltage  and  the  current 
under  our  specific  conditions,  an  order  of  magnitude  estimation 
has  been  conducted.  The  test  conditions  are  as  follows: 

gas  medium  =  nitrogen 

arc  length  =  1.4  cn  >  L  >  .6  cm  (cross  flow  configuration 
an  average  arc  length  of  1  cn  is  used  in  the 
calculation) 

gas  stream  static  pressure  =  p  ^  .07  atm  ~  53  mm  Eg 

gas  strean  static  temperature  =  T  ~  120 0°K 

gas  strean  Mach  number  =  M^,  ~  2.5 

The  literature  survey  reveals  that  in  Ref.  (29)  the  nost 
similar  case  to  our  test  conditions  was  investigated.  The  arc 
was  discharged  between  a  point  cathode  and  a  ring-shaped  anode 
in  coaxial  arrangement  and  were  35  sn  apart.  The  gas  medium  was 
argon  with  flow  Mach  number  3.  The  ambient  pressure  was  16  ran 
Eg.  A  minimum  arc  current  of  150  amps  and  a  maximum  arc 
voltage  of  about  60  volts  were  observed. 

PREDICTION  OF  MAINTENANCE  VOLTAGE 

In  order  to  extend  the  above  results  to  our  experiment, 
several  transformations  have  to  be  made  corresponding  to  our 
test  conditions. 

Transformation  of  the  Medium 

The  influence  of  the  flow  media  on  the  arc  discharge  has 
been  studied  by  John  and  Bade  (Ref.  30).  They  investigated 
different  flow  media  at  a  mass  flow  rate  of  about  0.5  gr/sec 
and  at  a  static  pressure  of  about  1  atm.  The  investigations, 
for  a  parallel  flow  configuration,  showed  that  the  smaller  the 
molecular  weight  of  the  medium  the  higher  the  voltage  required 
to  maintain  the  arc.  In  the  case  of  nitrogen  and  argon,  the 
required  voltage  ratio  is  about  3  to  2  (nearly  independent  of 
the  arc  current).  The  maximum  arc  voltage  given  in  Ref.  (29) 
should  be  corrected  accordingly  to  a  nitrogen  stream  which 
gives  V  %  90  volts. 
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Transformation -of  the  Mass  Flow  Rate 

For  air  at  a  static  pressure  of  about  1  atm,  the  influence 
of  the  mass  flow  rate  has  also  been  studied  by  John  and  Bade 
(Ref.  30).  From  the  observation  of  a  cross  flow  configuration, 
the  necessary  voltage  to  maintain  the  arc  is  increased  if  the 
mass  flow  rate  is  increased.  More  specifically,  the  voltage 
is  increased  by  about  20%  (independent  of  the  arc  current)  if 
the  mass  flow  rate  is  doubled. 

In  our  case,  a  rough  estimation  of  the  mass  flow  rate 
passing  through  the  arc  gives  m  <  3.2  gr/sec,  where  the  arc 
radius  has  been  assumed  to  be  0.T  cm.  This  is  about  1/50  of 
our  total  mass  flow  rate  of  the  nozzle.  In  view  of  this,  the 
arc  voltage  should  be  corrected  to  V  £  200  volts. 

Transformation  of  the  Arc  Length 

Busz  and  Finkelburg  (Ref.  31)  investigated  the  effect  of 
the  electrode  gap  for  an  argon  arc  at  200  amps.  Their  results 
showed  that  the  larger  the  arc  length  the  larger  the  required 
voltage  to  maintain  the  arc.  If  it  is  assumed  that  the  same 
qualitative  behavior  exists  for  a  nitrogen  arc,  then,  according 
to  Busz  and  Finkelburg  (Ref.  31)  the  arc  voltages  compare  like 


^3  5mm 
^6ram 


&  2 


Hence,  a  corrected  value  for  the  arc  voltage  should  be  about 
100  volts. 

Transformation  of  the  Ambient  Pressure 

In  Ref .  (32) ,  the  relation  between  arc  voltage  and  ambient 
pressure  is  given  in  detail.  It  can  be  observed  that  the  arc 
voltage  decreases  as  the  pressure  decreases.  However,  it  may 
be  found  that  this  decrease  of  the  arc  voltage  for  low  pressure 
arcs  is  less  substantial  than  in  case  of  high  pressure  levels. 
In  view  of  Ref.  (32),  it  may  be  assumed  that  no  transformation 
is  needed. 

The  maintenance  voltage  of  the  arc  in  our  case  may  then  be 
said  to  be  about  100  volts.  It  is  known  from  experiments  (Ref. 
32)  that  in  order  to  establish  a  stable  arc  it  is  required  that 
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V  £  2V 
o  a 


Therefore,  the  required  power  source  voltage,  V0,  is  about  200 
volts  in  our  case,  while  the  arc  power  supply  should  be  15KW 
approximately.  Using  the  voltage  balance  of  the  circuit,  the 
outer  circuit  resistance  may  be  calculated  to  be  about  0.7  ohms. 

ESTIMATION  OF  BREAKDOWN  VOLTAGE  RANGE 

Another  property  which  has  to  be  estimated  is  the  breakdown 
voltage.  The  estimation  is  based  on  the  results  given  in  the 
classical  work  of  Cobine  (Ref.  18)  for  the  following  test 
parameters: 

electrode  gap:  0.6  cm 

static  pressure:  53  mm  Hg 

Figure  IV-1  (reproduced  from  Ref.  (18))  may  be  used  to  estimate 
the  breakdown  for  the  arc  to  23,000  volts  approximately. 

As  it  is  pointed  out  by.  Cobine  (Ref .  18) ,  this  breakdown 
voltage  may  be  decreased  considerably  if  there  is  an  insulator 
connecting  the  electrodes.  This  is  due  to  the  phenomenon  of 
collection  of  surface  charge.  In  addition,  the  presence  of 
any  gas  layer,  moisture  or  dust  again  may  decrease  the  break¬ 
down  voltage.  Under  these  conditions  it  is  recommended  (Ref.  18) 
to  allow  a  factor  of  3  in  the  gap  length.  Moreover,  Howell  (see 
the  same  reference) ,  has  found  that  at  high  pressures  the 
breakdown  voltage  can  be  reduced  by  a  factor  of  4  if  the  surfaces 
of  electrodes  are  rough.  Therefore,  the  minimum  breakdown 
voltage,  assuming  rough  electrodes  and  surface  charge,  could  be 
as  low  as  approximately  400  volts. 


Breakdown  Voltage,  kilovolts 


Figure  IV-1.  Breakdown  Voltage  between  Cylindrical 
Rods  at  Atmospheric  Pressure 


APPENDIX  V 

ANALYTICAL  DETERMINATION  OF  THE  ARC  CHARACTERISTIC 


"Employing  the  Ellehbaas  Energy  Equation,  an  approximate 
method  to  determine  the  operation  parameter  of  an  electric  arc 
discharge  hhs  been  established  by  Msecker  (Ref.  19)  using  the 
Schmitz’s  approximation.  The  analysis  given  herein  connects 
the  characteristics  of  the  electric  discharge  with  the  parameter 
of  the  outer  circuit.  It  should  be  noted  that  Maecker’s 
original  analysis  was  derived  for  a  contained  cylindrical  arc. 
Since  no  theory  exists  in  predicting  the  arc  characteristic 
corresponding  to  our  case,  it  is  felt  that  Maecker’s  theory 
may  be  used  in  order  to  obtain  at  least  some  qualitative 
information. 

The  voltage  balance  for  the  whole  circuit  is  (Fig.  V-l): 


VQ  -  RaI  -  RCI  =  0  (1) 


where 

VQ  «  battery  or  power  source  voltage 

R  =  resistance  of  the  arc 
a 

R  =  resistance  of  the  outer  circuit 
c 

I  =  current  of  the  circuit 

Under  the  assumption  of  a  cylindrical  arc  the  electric  field 
strength  may  be  written  as 


IS 


where  1/  stands  for  the  arc  length.  Then,  Equation  (1)  may  be 
rewritten  as: 


90 


at  point  of  V-L<w 

operation 

where  Ea(I)  is  the  arc  characteristic,  or  the  current-voltage 
relationship  of  the  arc. 

Neglecting  any  mass  flow  as  well  as  radiation,  the  Ellenbaas 
Energy  Equation  in  cylindrical  coordinates  takes  the  following 
form: 


a(T)Ef  +  -  — [rk(T)~ ]  =  0  (2) 

r  dr  dr 


where 

r  =  arc  radius 

o(T)  =  electrical  conductivity 
k  =  thermal  conductivity 
.T  =  temperature  of  the  gas 

Introducing  the  heat  conduction  function, 

T 

S  =  /  k(T)dT  (3) 

O' 


then,  the  energy  equation  takes  the  form 


c(T)E^  +  I  A[r  £S]  =  o  (4) 

a  r  dr  dr 


Yn  R  I  V 

^  _a  _  $  (i) 

V  V  V  a 
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In  general,  the  electrical  conductivity  oCO  varies  rather 
significantly  with  gas  temperature,  hence  the  material  function 
is  a  =  <r(S)  -  Usually  this  is  determined  experimentally  and 
given  in  graphs.  In  view  of  the  substantial  temperature  differences 
between  the  core  of  the  arc  discharge  and  the  boundary,  the 
material  function  may  be  approximated  by  two  linear  segments 
dividing  the  arc  into  two  regions,  (Fig*  V-2).  Let  ex*  be 
the  electrical  conductivity  corresponding  to  a  point  of  the 
approximating  line  segment  (Fig.  V-2)  then  o  =  c(S)  may  be 
expressed  as: 


o*  =  0  for  0  <  S  <  S 

a 

a*  =  a(S  -  S  )  for  S  <  S  <  S 
a  a  o 


where  a  is  the  slope  of  one  of  the  line  segments.  As  mentioned, 
the  above  approximation  divides  the  arc  into  two  regions,  namely, 
the  outer  region  without  electrical  conductivity  and  the  inner 
region  with  electrical  conductivity.  These  two  regions  can  be 
considered  separately  as  follows: 

The  Region  with  Electric  Conductivity 

In  this  region  a*  =  a(S  -  S  ) ,  defining  x  =  rE  v/a~,  we 
find  that  Equation  (4)  may  be  deduced  to  a 


a*  +  —  a*  +  a*  =  0 

XX  x  X 


(5) 


where  subscript  x  denotes  the  differentiation  with  respect  to 
x.  This  equation  is  a  zeroth  order  Eessel  Equation.  Since 
o*(x)  has  to  be  finite  for  all  x  and  J  (x  =  0)  =  1,  we  have  - 


a*(x)  =  o*  JQ(x) 


(6) 


The  last  term  follows  simply  from  the  geometry  of  Fig.  V-2.  We 
note  that  for  x  =  xQ  =  2.405  the  Bessel  function  reaches  its 
first  zero;  in  that  case  o*(x  =  2.405)  =  0  and  S  =  S  . 

a 
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The  Region  Without  Electric  Conductivity 

In  thi &  region  the  energy  equation,  i.e.,  Equation  (4) 
becomes  C.j=  x  dS/dx.  In  order  to  determine  the  integration 
constant  C  ,  we  require  that  the  slopes  dS/dx  for  the  outer  and 
inner  regions  to  be  matched  at  the  intersection.  From  the 
solution,  of  the  inner  region,  we  have 


dS! 

dx 


o 


<So  -  Sa> 


d(J  rx]) 


dx 


-  -  <so  -  WV 


Multiplying  this  by  xQ  =  2.405,  we  have 

»» 


x 


o 


.dS 

dx 


=  -  1.248(Sq  -  Sa) 
o 


Matching  this  result  with  the  outer  region  solution,  we  obtain 


C  =  x  —  =  -  1.248(S  -  S'  ) 

1  dx  o  a 


Now  a  second  integration  of  dS  =  -  1.248(SQ  -  Sa) —  gives,  using 
the  boundary  condition  that  for  x  =  xQ,  S  =  S&  x  . 


1.248  In.— 


x 


o 


(7) 


Replacing  x  by  its  definition  x  =  rE ^/oT,  and  requiring  that 
r  =  r  and  S  =  0  at  the  wall  the  following  relation  is  obtained: 


Jk 

1.248  In  I— - — (EyE" 
2.405  a 


(8) 
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where  r  is  the  arc  tube  radius  in  Maecker's  original  text. 
The  electrical  field  strength  is 


E  =  2.405 
a 


_i_  .  er.248(so  -  sa) 

V5" 


O) 


Writing  Equation  (7)  for  x  =  r  and  Xq=  r  ,  the  following  equation 
is  obtained:  a 


Pa  =  (2ir)  (1.248)  (SQ  -  S&) 

=  7.84(S0  -  Sa)  (13) 

Oh  the  other  hand,  Equation  (8)  after  some  manipulation,  can  be 
rewritten  as 

E  =  •  exp.[ - - ]  (14) 

raV/cT  1.248(So  -  S&) 
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so  that, 


P  0*  S 

—  =  jl  =  7.27  r  —  exp .  [ - 2 - ]  (15) 

Ea  aV5"  1.248(SG  -  Sa) 


It  can  be  observed  that  for  constant  arc  power  and  a  given  gas, 
the  values  of  SQ  and  Tc  are  independent  of  the  arc  tube  radius. 

On  the  other  hand,  it  is  evident  that  the  arc  current  varies 
proportionally  and  the  electric  field  varies  inverse  proportionally 
with  the  arc  radius. 

To  approximate  the  material  function  0(S)  with  the  line 
segment  of  Fig,  Y-2,  two  filling  factors  f  and  fx  were  defined 
by  Maecker  (Ref.  19)  such  that: 


1  *° 

s-i-  /  crdS  ; 


a  S  „ 
o  o  0 


f  =  - . 

1  S0/°odS  o 
0 


J  0SdS 


(16) 


Now  the  Equations  (13) ,  (14)  and  (15)  may  be  rewritten  to  obtain: 


Pa  =  6irzf  x  (T)Sq(T)  ,  [kw/cm]  (17) 


3  94  / (T) 

E  =  f  (T)  / - 2 -  exp.  [ - ± - ] ,  [volt/cm]  (18) 

a  r  1  V  2f  (T),.0„(T)  ,3zf  (T) 

a.  »  O  2 

I  =  7.27  rJ2f(T)S0(T)oo(T)  exp  .  [ - ,  [amp]  (19) 

V  3zf 


where  z  =  1.248.  These  equations  enable  us  to  find  the  arc 
characteristics  provided  S0(T),  0q(T),  f (T)  and  fr(T)  are 
known  as  functions  of  the  arc  center  temperature.  In  addition 
it  is  also  necessary  to  know  the  arc  tube  radius. 
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To  demonstrate  the  method  of  determining  the  arc  character¬ 
istics  from  measured  values  of  arc  current  and  arc  voltage,  the 
following  test  data  will  be  considered: 


Arc  Current 

Arc  Voltage 

Arc  Power 

I 

V 

P 

a 

a 

[amps] 

[volts] 

[kw] 

26 

115 

2.99 

64 

25 

1.60 

104 

57.5 

5.98 

94 

60 

5.64 

102 

47.5 

4.84 

148 

40 

5.92 

78 

0 

0 

These  data  are  given  in  Fig.  V-3  assuming  that  the  average  arc 
length  is  1  cm.  As  it  has  been  mentioned  the  arc  characteristic 
may  be  determined  with  the  aid  of  the  Equations  (17),  (18)  and 
(19).  The  unknown  properties  S(TC),  o(T?) ,  f  (T?)  and  f(Tc)  are 
given  in  Ref.  (19)  for  nitrogen  as  functions  5f  the  axis  temper¬ 
ature  of  the  arc.  The  temperature  at  the  arc  axis  are  not  known 
a  priori.  However,  they  may  be  estimated  as  follows. 

Assuming  a  linear  decrease  of  the  arc  temperature  from  its 
axis  value  Tc  to  its  edge  value  Ta  (here  the  edge  is  the  boundary 
between  the  arc  and  the  external  flow  field) .  The  average 
value  for  the  arc  temperature  may  be  defined  as 


T  = 


T  +  T 
c  a 


(20) 


So  that  the  axis  temperature  may  be  written  as 


T  =  2T  -  T 
c  a 


(21) 
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The  average  plasma  temperatures  are  given  as  functions  of  the 
arc  power  in  Ref. (33).  It  should  be  noted  that  even  though 
the  temperature  values  given  in  Ref.  (33)  are  corresponding  to 
an  argon  plasma,  the  graphs  as  shown  in  Fig.  V-4  should  still 
give  qualitative  representation  of  the  temperature,  at  least 
within  the  limits  of  our  approximation.  The  effect  of  the  mass 
flow  rate  upon  the  average  temperature  may  be  neglected  (Ref .  34) . 

Now  the  edge  temperature  of  the  arc  has  been  assumed  to  be 
the  static  temperature  of  the  free  stream  (T  ft  1200°K) .  It 
should  be  noted  however,  that  this  temperature  may  be  somewhat- 
low  as  compared  to  the  actual  case.  In  Fig.  V-3  the  calculated 
power  line  corresponding  to  that  edge  temperature  of  the  arc 
is  shown  with  the  measured  data.  The  agreement  is  fair .  The 
discrepancy  can  be  explained  simpiy  by  the  assumption  made 
concerning  the  temperature  at  the  arc  boundary  being  the  static 
temperature  of  the  free  stream.  Obviously  the  assumed  value 
was  too  low;  this  possibility  has  already  been  mentioned.  In 
view  of  the  Equations  (17)  and  (21),  it  might  be  observed  that 
it  is  possible  -  using  an  iterative  procedure  -  to  determine 
the  correct  value  of  the  edge  temperatures.  Fig.  V-4  shows  ■ 
the  calculated  arc  characteristics  for  different  arc  radii  as 
•well  as  the  measured  data.  Again,  the  correct  arc  radius  may 
be  obtained  by  an  iterative  procedure. 

It  should  be  pointed  out  here  that  Maecker's  method  was 
derived  for  a  contained  cylindrical  arc,  hence,  the  results 
cannot  be  safely  claimed  to  be  directly  applicable  to  our 
case  in  a  strict  manner.  However,  since  we  are  interested 
mainly  in  obtaining  some  qualitative  data  concerning  arc 
center  and  arc  boundary  temperature  as  well  as  the  arc  radius, 
the  way  of  obtaining  these  quantities  as  stated  above  seems 
justifiable.  .  One  may  observe  that  the  measured  arc  character¬ 
istic  is  most  closely  represented  by  the  arc  characteristic 
of  an  arc  operating  in  a  tube  of  approximately  0.1  cm. 
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APPENDIX  VI 


DIFFICULTIES  OF  THE  EXPERIMENT 


Information  on  pressure  distribution  in  a  flow  field  dis¬ 
turbed  by  electric  discharge  is  not  available  in  the  literature. 
The  magnitudes  of  the  pressures  and  their  affected  areas  could 
only  be  evaluated  by  using  the  simplified  theories  established 
for  thrust  vector  control  by  secondary  injections.  Different 
opinions  in  theoretical  considerations,  as  discussed  in  Appendix 
I,  made  it  necessary  to  conduct  an  experimental  investigation. 
However,  the  mechanism  of  creating  the  flow  field  disturbance 
by  an  electric  arc  discharge  was  not  well-understood;  for 
example,  the  breakdown  voltage  of  an  electric  arc  could  not 
be  analytically  determined  in  a  supersonic  flow  field  (Appendix 
IV) .  The  additional  material  problems  for  both  the  nozzle 
section  and  the  electrode  insulations  also  needed  to  be  explored 
In  order  to  solve  all  these  major  problems,  a  water-cooled  axi- 
symmetric  copper  nozzle  was  built  (Fig. 10)  for  preliminary 
investigation  in  order  to  determine  the  design  criteria. 

PRELIMINARY  TEST 

The  first  test  was  successfully  completed  with  the 
following  information  confirmed: 

1)  The  water-cooled  copper  nozzle  designed  at  a  marginal 
strength  by  considering  steady  state  heat  transfer 

of  conduction  and  convection  sustained  the  high 
temperature  environment  for  ten  to  thirty  seconds 
without  problem . 

2)  The  arc  can  be  established  between  a  pair  of  cross-flow 
electrodes  of  0.25  inch  apart  at  a  breakdown  voltage  of 
375  volts  for  a  nitrogen  stream  of  Mach  number  2.8 
with  static  pressure  of  2  psia  and  stagnation  tem¬ 
perature  of  4000°R. 

3)  The  electrode  configurations  affected  the  breakdown 
voltage  as  well  as  the  effective  area.  Since  the 
highest  breakdown  voltage  was  considered  within  the 
practical  range,  it  was  felt  that  the  largest  affected 
area  should  be  more  useful  for  thrust  vector  control 
applications.  The  largest  effective  area  was  obtained 
for  the  configuration  that  the  two  electrodes  were 
perpendicular  to  the  flow  direction. 
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4)  An  electric  arc  can  cause  the  static  pressure  rise  in 
a  supersonic  flow  field.  The  arc  can  be  maintained 
without'  the  confinement  of  a  magnetic  field. 

The  difficulties  encountered  in  the  test  for  the  preliminary 
investigations  were  as  follows: 

Flow  Separations 

The  nozzle  was  designed  for  a  back  pressure  of  1  psia 
which  was  originally  expected  from  the  test  facility.  Leakage 
problems  in  the  vacuum  system  were  not  able  to  be  solved  during 
the  time  allocated  for  this  test.  Thus  the  pressure  distributions 
(Figs.  26,  27  and  28)  at  the  nozzle  exit  should  not  be  considered 
relevant  to  the  purpose  of  this  test.  A  back  pressure  of  2  psia 
must  be  used  as  the  design  parameter  of  the  future  test  nozzles. 

Material  Problems 

— ■  ■  — — —  m  i  ■in———— 

The  insulation  material  between  the  electrodes  and  the  copper 
was  severely  damaged  as  shown  in  Fig.  29.  Micarta,  silicon 
phenolic,  and  boron  nitride  were  used.  Micarta  was  melting  and 
running  with  the  flow  almost  instantly.  Silicon  phenolic  became 
a  good  electric  conductor  at  high  temperature.  Boron  nitride 
was  ablated  at  a  comparatively  slow  rate  and  was  considered 
the  best  insulation  material  for  this  purpose. 

The  method  of  electrode  protrusion  for  protecting  the 
insulation  was  tested  by  extending  the  electrodes  approximately 
one-half  of  the  boundary  layer  displacement  thickness  into  the 
free  stream.  The  insulation  was  comparatively  better  protected 
but  the  flow  field  was  also  disturbed  substantially.  Fig.  IV-1 
shows  the  pressure  distributions  caused  by  electrode  protrusion 
with  and  without  discharge.  It  was  thus  felt  that  the  electiode 
protrusion  did  not  have  any  advantage  and  should  not  be  used 
for  further  testing. 

Pressure  Response  Rate 

Transducers  and  manometers  were  used  for  pressure  measure¬ 
ments.  Due  to  the  distance  between  the  apparatus  and  the  instru¬ 
mentation  panel,  a  finite  time  period  was  required  to  reach  a 
steady  state  reading.  The  pressure  response  history  was' 
recorded  and  is  shown  in  Fig.  VI-1  for  both  transducers  and 
manometers.  It  can  be  seen  that  the  transducers  required  0.8 
seconds  and  the  manometers  required  3.0  seconds  to  reach  the 
true  pressure  values  of  the  flow  field. 


Data  Uncertainty 


The  material  problem  required  the  discharge  time  to  be 
as  short  as  possible  in  order  to  produce  minimum  damage  which 
influenced  the  flow  field  pressure  readings.  On  the  other 
hand,  the  pressure  response  rate  requires  longer  time  in  order 
to  register  the  steady  state  readings.  Since  these  two  demands 
could  not  be  satisfied  simultaneously,  a  compromised  time 
interval  had  to  be  selected.  Three  seconds  fox-  electric 
discharge  was  decided  for  the  future  test  with  pressure  data 
taken  after  the  discharge  for  comparing  with  that  before  dis¬ 
charge.  If  any  disagreement  occurred  between  the  after  and 
before  discharge  values,  the  pressure  differences  obtained  by 
subtracting  the  after  discharge  value  from  the  during  discharge 
value  should  be  used  in  order  to  give  a  more  conservative  result 
for  side  force  evaluation. 

Practice  of  Side  Force  Evaluation 


frith  the  static  pressure  data  from  the  preliminary  test, 
it  was  possible  to  evaluate  the  generated  side  force  for  both 
practice  and  curiosity.  The  practice  was  to  determine  the 
actual  difficulty  which  might  be  involved.  The  curiosity  was 
to  check  the  different  opinions  of  Appendix  I ,  since  an  order 
of  magnitude  difference  of  the  generated  side  force  was  in 
question. 

It  was  felt  that  the  pressure  measurements  were  not  as 
comprehensive  as  needed  in  order  to  give  a  quantitative  side 
force  evaluation.  However,  since  the  preliminary  investigation 
had  to  provide  the  various  combinations  for  electrode  configura¬ 
tions,  the  limited  pressure  data  were  used  for  the  order  of 
magnitude  comparison  to  verify  the  different  opinions  of  the 
analytical  considerations . 

The  procedure  of  Appendix  VII  was  used  for  the  side  force 
evaluation.  With  very  limited  numbers  of  data  points  on 
pressure  distributions,  only  the  order  of  magnitude  of  the 
side  force  was  able  to  be  evaluated  and  agreed  well  in  the  order 
of  magnitude  with  Broadwell's  modified  linearized  theory 
(Appendix  I).  This  information  was  exciting  because  if  could 
bring  the  thrust  vector  control  by  electric  discharge  comparable 
with  the  secondary  injection  scheme  even  with  present  technology 
on  power  sources,  namely  chemical  batteries.  The  determining 
factor  of  the  success  of  this  new  concept,  however,  lies  on 
a  detailed  verification  of  the  pressure  measurements  in  the 
entire  flow  field  especially  in  the  vicinity  downstream  of  the 
electrodes . 


DETAILED  EXPERIMENT 

From  the  experience  of  the  preliminary  test,  a  rectangular 
nozzle  was  designed  in  order  to  satisfy  the  following  purposes: 

1)  Detailed  pressure  measurements 

2)  Visualization  of  the  flow  field  caused  by  the  electric 
discharge 

3)  Hot  and  cold  boundary  layer  effects  at  the  discharge 

4)  Electrode  configuration ‘effect  on  the  high  pressure 
areas . 

The  nozzle  consisted  of  three  parts  as  shown  in  Fig.  18.  The 
bottom  plate  could  be  replaced  easily  for  boundary  layer  and 
electrode  configuration  investigations.  Three  different  materials 
were  used  for  the  bottom  plate: 

Boron  Nitride 


From  the  insulation  material  experience,  the  boron  nitride 
appeared  to  be  the  most  attractive  choice  because  of  its  high 
thermal  shock  resistance.  A  plate  was  built  with  heavy  concen¬ 
tration  of  pressure  taps  in  the  vicinity  of  the  electrodes. 
Pressure  checks  indicated  very  serious  leakage  problems  for 
all  pressure  taps.  Several  methods  were  tried  for  preventing 
the  leakage,  only  the  metallic  coated  boron  nitride  plate  appeared 
to  be  possible.  However,  the  cost  and  time  needed  for  this 
process  made  it  impractical  because  the  number  of  pressure  taps 
and  different  configurations  required. 

Silicon  Phenolic 


Silicon  phenolic  appeared  to  be  the  second  to  the  best  choice, 
if  the  electrodes  were  insulated  with  boron  nitride.  This  plate 
was  made  and  tested.  The  flow  field  disturbance  by  electric 
discharge  was  clearly  observed  through  quartz  windows  downstream 
of  the  electrodes.  The  only  problem  was  the  deformation  of  the 
plate  at  high  temperature.  Severe  ablation  of  the  phenolic 
material  occurred  in  the  vicinity  of  the  electrodes.  The  static 
pressure  after  discharge  did  not  return  to  the  values  before 
discharge  as  shown  in  Fig.  33.  Nevertheless,  the  pressure  dis¬ 
tribution  is  different  from  the  interpolated  results  of  the 
axisyrametric  copper  nozzle  where  the  sharp  pressure  drop  right 
behind  the  pressure  peak  could  not  be  determined.  Side  force  .. 
evaluation  was  made  by  using  the  pressure  difference  during  and 
after  discharge.  The  evaluated  side  force,  however,  was  an 
order  of  magnitude  smaller  than  that  obtained  from  the  preliminary 
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test  and  agreed  well  with  the  theoretical  results  of  Section  II. 
It  was  felt  that  the  only  way  to  confirm  the  magnitude  of  the 
side  force  should  be  achieved  from  the  test  results  where  the 
pressure  measurements  after  discharge  returned  to  that  of 
before  discharge. 

Water-cooled  Copper 

Water-cooled  copper  appeared  to  be  free  from  deformation 
and  was  used  for  the  bottom  plate  with  boron  nitride  as  electrode 
insulators.  The  discharge  was  established  from  the  cathode  to 
copper  then  to  the  anode.  This  process  transformed  a  large 
quantity  of  the  supplied  electric  energy  to  the  cooling  water. 
Ceramic  spray  of  the  copper  plate  was  used  and  proved  to  be  the 
best  method  for  this  purpose.  Pressure  distributions  were 
measured  as  shown  in  Fig.  36.  The  pressure  readings  after  and 
before  discharge  were  identical,  but  the  area  with  pressure  rise 
during  discharge  was  small.  Due  to  the  need  of  the  boron  nitride 
insulators,  no  pressure  tap  was  available  between  the  electrodes. 
Nevertheless,  it  is  obvious,  the  peak  pressure  rise  dropped 
sharply  along  the. flow  direction.  The  side  force  evaluation 
confirmed  the  theory  of  Section  II. 

Silicon  Phenolic  Axisymmetric  Nozzle 

A  silicon  phenolic  axisymmetric  nozzle  was  also  tested  for 
the  purpose  of  confirming  the  rectangular  nozzle  results  and 
comparing  the  cold  boundary  layer  effect.  The  pressure  distri¬ 
bution  is  shown  in  Fig.  39.  The  evaluated  side  force  again 
confirmed  the  theoretical  results  of  Section  II. 

CONCLUDING  REMARKS 

It  was  felt  that  an  additional  test  should  be  conducted 
by  using  the  rectangular  nozzle  with  ceramic  sprayed  copper  plate 
However,  more  pressure  measurements  with  faster  responses  would 
be  desired.  A  modification  was  needed  by  re-fabricating  a  new 
copper  plate  and  re-installing  the  pressure  transducers  closer 
to  the  test  nozzle. 

The  project  was  terminated  for  the  following  reasons: 

1)  Additional  time  and  funds  would  be  required  for  further 
testing  which  appeared  to  conflict  with  the  higher 
priority  projects  scheduled  for  the  test  facility  in 
the  near  future. 


2)  The  material  problem  is  very  severe  from  the  system 
applications  point  of  view.  This  could  not  be  solved 
even  with  more  precise  test  data  available. 
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APPENDIX  YU¬ 


AN  EXAMPLE  OF  THE  SIDE  FORCE  CALCULATION 


An  example  to  illustrate  the  side  force  evaluation  is 
presented  here  for  the  purpose  of  indicating  the  attainable 
accuracy  in  the  calculation.  The  test  results  obtained  from 
the  axisymmetric  silicon  phenolic  nozzle  are  used  for  the 
numerical  calculation.  The  pressure  tap  distributions  onto 
the  projected  area  are  shown  in  Fig.  VII-1.  The  x-axis  is  along 
the  flow  direction  and  the  y-axis  normal  to  the  flow  direction. 
Since  the  electrodes  are  symmetric  to  the  nozzle  centerline,  the 
pressure  distributions  are  also  assumed  to  be  symmetric  to  the 
same  centerline. 

Pressure  readings  before,  during  and  after  discharge  have 
been  all  recorded.  Fig.  VII-2  shows  the  reading  along  nozzle 
centerline  (yQ-seetion) .  Figure  VI 1-3  and  Figure  VI 1-4  show  Ifae 
readings  along  x^section  and  x2 -section,  respectively. 

Similarly,  readings  in  sections  x3,  x4,  x5,  y  ,  and  y2  can  be 
obtained.  Due  to  the  material  deformation  under  high  temperature 
environment,  it  is  noted  that  the  pressure  readings  before  and 
after  discharges  are  not  consistent.  In  order  to  be  able  to 
study  the  feasibility  of  this  concept,  the  pressure  rise  caused 
by  electric  discharge  is  assumed  to  be  the  pressure  difference 
between, during  (P^urin^  and  after  (Pafter)  discharges.  The 
side  force  is  then  evaluated  as 


Fs  =  2  /  «Wing  -  Pafter)dA 
A 


where  A  is  one-half  of  the  affected  area  or  the  area  behind  the 
shock.  Dividing  the  affected  area  A  into  some  reasonably  small 
subdivisions,  Mi,  and  since  these  subdivisions  are  sufficiently 
small,  we  can  assume  that  the  pressure  change  is  constant  within 
each  subdivision.  In  some  high  pressure  gradient  regions  a 
finer  partition  is  used.  The  mean  value  of  the  pressures  is 
chosen  at  each  subdivision.  Those  mean  values  are  obtained  from 
Figures  VII-2,  VII-3,  and  VII-4.  Let, 


M^  =  the  area  of  each  subdivision, 


then 


(Pduring  "  Pafter)iMi  =  the  side  force  Produced  in  ^i* 


Summing  up  all  the  affected  areas,  Mi,  we  obtain 


F  =  22 
s 

i 


(P 


daring 


^after^i^i  =  si£ie  force. 


The  procedure  of  the  side  force  calculation  for  the  silicon 
phenolic  nozzle  is  as  follows: 

1  5 

1.  Choose  Mi  =  —  in  as  shown  in  Fig.  VTI-5. 

64 

2.  Find  all  the  pressure  rises 


^i  ^during  ~  Pafter^i 


3.  Trace  out  the  isobaric  area.  Namely,  to  determine 
the  number  of  small  blocks,  Ni,  with  same  values  of 
pressure  rise  APi- 

4.  Then,  we  can  write 


F„  =  2  2  (P 


during  "  ^after^^i 


=  22  APi(NiMi) 
i 


=  2(Mi)  2  (AP^CNj.) 

i 

which  gives 


F  =  2  x  (— )[(0.1)(18.5)  +  (0.2)  (15.5) 
s  64 

+  (0.25) (3)  +  (0.3) (9)  +  (0.8) (0, 5) 

+  (-  0.2) (1)  +  (-  0.5) (1.5)  +  (-  0. 8) (0. 5) ] 

i 

=0.26  lbs. 
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Figure  VTI-1.  The  Axisysimetric  Silicon  Phenolic  Nozzle  with  0.5  in. 
Electrode  Distance  in  Cross-Flow 


Ap,  psi  Static  Pressure,  psia 


Figure  VI 1-2.  Pressure  Readings  along  Center-Line  of  the 
Axisymmetric  Silicon  Phenolic  Nozzle. 


□0 


Ap,  psi  Static  Pressure,  psia 


Figure  VII-3,  Pressure  Readings  along  X, -Section 

(y  =  0  is  the  projection  of  Y  -Section) 


O :  Before  Discharge 
A :  During  Discharge 
□  :  After  Discharge 
AP  «  Prflir,ntf  -  P 


during  “  ^after 


NOTE:  Those  indicated  values 
changes  in  that  area 


Figure  VII-5.  pressure  Change 
Affected  Area 
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